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1. Introduction 
Consider a polynomial line field on the real (a;, y)-planc 

(1) ^ = ^ l [ X,V \ , P n ,Qn ~ polynomials ,degP„,Q„ < n. 

dx Q n [x,y) 



H(n) = uniform bound for the number of limit cycle of (1) 



One way to formulate the Hilbert 16-th Problem is the following: 

Hilbert 16-th Problem (HP). Find an estimate for H{n) for any n 6 Z+. 

We shall discuss problems related to the following: 

Existential Hilbert 16-th Problem (EHP). Prove that H{n) < oo for any 
n e Z+. 

The problem about finiteness of number of limit cycles for an individual polyno- 
mial line field ([l]) is called Dulac problem since the pioneering work of Dulac who 
claimed in 1923 to solve this problem, but an error was found by Ilyashenko. 

The Dulac problem was solved by two independent and rather different proofs 
given almost simultaneously by Ilyashenko B and Ecalle jEj. However, both proofs 
do not allow any generalization to solve Existential Hilbert Problem. 

Consider the equation (Q) for different polynomials [P n {x, y), Q n (x, y)) as the 
family of line fields on R 2 depending on parameters of the polynomials. Using a 
central projection ir : § 2 — > R 2 and homogcnuity with respect to parameters of the 
equation ([jj) (line fields XP n (x,y)/XQ n (x,y) and P n (x,y)/Q n {x,y) for any A ^ 
are the same) one can construct a finite parameter family of analy tic line fields 
on the shpere S 2 with a compact parameter base B (see e.g. |jY2[ for details). 
After this reduction Existential Hilbert Problem becomes a particular case of the 
following 

Global Finiteness Conjecture (GFC). (see e.g. [Qj For any family of line 
fields on § 2 with a compact parameter base B the number of limit cycles is uniformly 
bounded over all parameter values. 



We refer the reader to the volumes || and IY2] where various development of 



these and related problems are discussed. Families of analytic fields are extremely 



difficult to analyze. In the middle of 80's Arnold AAI | proposed to consider generic 
families of smooth vector fields on S 2 . A smooth analog of Global Finiteness Con- 
jecture is the following 
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Hilbert- Arnold Problem (HAP). (e.g- JTYSfl j Prove that in a generic finite 
parameter of vector fields on the sphere S 2 with compact base B , the number of 
limit cycles is uniformly bounded. 

Assume for a moment that a polynomial (or a generic smooth) vector field on 
the sphere § 2 has an infinite number of limit cycles. By the Poincare-Bendixon 
Theorem, any limit cycle should surround an equilibrium point and, since our vec- 
tor field has at most finitely many equilibria, there should be an infinite "nested" 
sequence around one of equilibria. Then those "nested" limit cycles have to accu- 
mulate (in the sense of Hausdorff metric) to a certain contour (polygon) consisting 
of equilibria (as vertices) and separatric curves (sides of that polygon) connecting 
them. Such objects are called polycycles. It turns out that a possible solution to 
Hilbert- Arnold Problem reduces to investigation of bifurcation of polycycles. Let 
us give several definitions. 

Definition 1. A poly cycle 7 of a vector field on the sphere § 2 is a cyclically ordered 
collection of equilibrium points p\ , . . . , pk ( with possible repetitions ) and different 
arcs 71, . . . ,7^ (integral curves of the vector field) connecting them in the specific 
order: the j-th arc jj connects pj with Pj+i for j = 1, . . . , k. 

Definition 2. Let {x = v(x, e)} ee B"j x G § 2 , be an n-parameter family of vector 
fields on § 2 having a poly cycle 7 for the critical parameter value e* . The poly cycle 
7 has cyclicity fi in the family {v(x, e)} e£ s n if there exist neighborhoods U and V 
such that S 2 D U D 7, B D V G e* and for any e G V the field w(-, e) has no more 
than fj, limit cycles inside U and fi is the minimal number with this property. 

Examples 1) In a generic n-parameter family, the maximal multiplicity of a 
degenerate limit cycle does not exceed n + 1, e.g. in codimension 1 a semistable 
limit cycle has cyclicity 2. Thus, the cyclicity of a trivial polycycle (a polycycle 
without singular points) in a generic n-parameter family does not exceed n + 1 . 

2) (Andronow-Leontovich, 1930s; Hopf, 1940s). A nontrivial polycycle of codi- 
mension 1 has cyclicity at most 1. 

3) (Ta kens, Bogdanov, Leontovich, Mourtada, Grozovskii, early 1970s-1993 (see 
Q, |KS| and references there)). A nontrivial polycycle of codimension 2 has cyclic- 



ity at most 2. 

Definition 3. The bifurcation number B{k) is the maximal cyclicity of a nontrivial 
polycycle occurring in a generic k-parameter family. 

The definition of B(k) does not depend on a choice of the base of the family, it 
depends only on the number k of parameters. 



Local Hilbert- Arnold Problem (LHAP) e.g.| [Yl | Prove that for any finite 
k, the bifurcation number B(k) is finite and find an upper estimate for B(k). 

It turns out that a solution to Local Hilbert- Arnold Problem implies a solution 
to Hilbert- Arnold Problem. 

Similarly to the generic smooth vector fields, in the case of analytic vector fields 
one can define so-called a limit periodic set JFP[ , [Q, |IY1 1 , which is either a poly- 



cycle or has an arc of equilibrium points^], and formulate 

Local Finiteness Conjecture (LFC) e.g.jRj] Prove that any limit periodic set 
occuring in an analytic family of vector fields on Jp has finite cyclicity in this family. 



generic vector fields can not have an arc of equilibrium points 



V. Kaloshin 



3 



Smooth vector fields are more flexible then analytic vector fields and easier to 
analyze. A strategy to attack Existential Hilbert Problem proposed by Arnold 
|AAI | (see also [ IK ) is first understand generic smooth vector fields and then try to 
apply developed methods to analytic vector fields. Let us summarize the discussion 
in the form of diagramm: 



EHP 



<= 



GFC 



<= 



LFC 



HAP 



<= 



LH AP 



FEC 



<= + 



RS & SP ? 



Figure 1. 



Now we shall formulate the Main Result of the paper. 

Definition 4. A singular (equilibrium) point of a vector field on the two-sphere is 
called elementary if at least one eigenvalue of its linear part is nonzero. A polycycle 
is called an elementary polycycle if all its singularities are elementary. 

The Local Hilbert- Arnold problem was solved under the additional assumption 
that a polycycle have elementary singularities only. 

Definition 5. The elementary bifurcation number E(k) is the maximal cyclicity of 
a nontrivial elementary polycycle occurring in a generic k-parameter family. 

From examples 2) and 3) above it follows that 

E(l) = l, E(2) = 2. 

Information about behavior of the function k i— > E(k) has been obtained recently. 
The First crucial step was done by Ilyashenko and Yakovenko: 

Finiteness Theorem (Ilyashenko and Yakovenko pY3j ]J For any n the elemen- 
tary bifurcation number E{n) is finite. 

Corollary 1. Under the assumption that families of vector fields have elementary 
singularities only the global Hilbert- Arnold conjecture is solved, i.e. any generic 
finite parameter family of vector fields on the sphere § 2 with a compact base and 
only elementary singularities has a uniform upper bound for the number of limit 
cycles. 

Main Theorem. For any k G Z + 

(2) E(k)< 2 25k \ 

This is the first known sufficiently general estimate for cyclicity of polycycle. The 
case of a polycycle consisting only one singular point with no arcs at all, is well 
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known. An elementary equilibrium point can generate limit cycles in its small 
neighborhood if it is a slow focus, that is the linearization matrix has a pair of two 



imaginary eigenvalues. This bifurcation was investigated by Takens [Ta 



Corollary 2. Under the assumption that all the polycycles are elementary the Main 
Theorem gives a solution to the Local Hilbert- Arnold problem. 

The Main Theorem is an improvement of Ilyashenko-Yakovcnko Finiteness The- 
orem. It is a great pleasure for the author to say that the paper of Ilyashenko- 



Yakovenko [IY3| was a corner stone for the present paper. In IY2] the authors 
made an extremely important step: they found a pass from bifurcation theory to 
singularity theory using the Khovanskii reduction method jKH] . We follow this pass 
up to some point and using some new ideas getting the first sufficiently general 
estimate for the cyclicity of polycycles. To make this paper readable we have to 



reproduce some points from [1Y3| and we are sorry for repetition, but we think that 



it is necessary for a better understanding. 

1.1. Three stages of the proof. The proof of the Main Theorem consists of three 



steps. Relation to the proof of the Finiteness Theorem [ [Y3 is discussed after this 
short description. 

Step 1. Normal forms for local families of vector fields and their integration In 
section || we use normal forms to establish an explicit form for the Poincare corre- 
spondence map near equilibrium points on the polycycle under consideration. In 
|IY3 it is shown that these maps satisfy Pfaffian (polynomial differential) equa- 



tions with coefficient of polynomials depending smoothly on the parameters of the 
family. As the result a basic system of equations for determination of limit cycles 
is obtained. 

Step 2. the Khovanskii reduction method In section ^ we discuss a variation 



of the Khovanskii method [Kh|. This method allows us to investigate systems 
of equations that involve functions satisfying Pfaffian equations. In section ^ we 
present a formal reduction from the basic system to a mixed functional- Pfaffian 



system which is done in [IY3] together with upper bounds for degrees of involved 
into the procedure polynomials. After application of the Khovanskii method to the 
mixed functional-Pfaffian system we obtain several chain maps, the maps of the 
form 

(3) x ~ (Pi, . . . , P„) o (x, f(x)J'(x), f M (x)), 

where P = (Pi, . . . , P n ) is a vector-polynomial given by its coordinate functions of 
known degree and / is a generic function. The problem of estimating the number of 
limit cycles reduces to estimating the number of regular preimages of some special 
points by the chain map. Special points form an open cone-like semialgebraic set 
K in the image. 

Denote by F the map F : x i— > P o (x, f(x), f'{x), . . . , f^ n '(x)) which is called 
the n-th jet of /. Denote by Lp the linearization of F at point x = 0. 

Step 3. Bezout's theorem for the Chain maps In section |5| we construct an 
algebraic set £ in the image of F (in the space of n-jets). If F is transversal to S, 
then the number of preimages of any point a from a set of special points K is the 
same for F and its linearization Lp at zero, namely, 

k- 

(4) #{x :PoF(x) = a} = #{a; : P o L F (x) = a} < JJdegP,-. 
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But if is a linear map and one can apply Bezout's theorem to estimate the right- 
hand side of the equality. This observation completes the proof of the Main Theo- 
rem. 

Let us discuss relation of this proof to the proof of the Finiteness Theorem by 
Ilyashenko & Yakovenko flY3|] . S tep 1 of both proof is the same. We just refer 
to appropriate statements in |IY3[ ] . Step 2 in this proof is slightly different for the 
one in ]IY3fl . After application of the Khovanskii method they obtain the same 
collection of chain maps of the form (||). However, they investigate the number of 
regular preimages of points in the image by the chain maps without any restriction 
on those points. In the present proof, using additional arguments in the Khovanskii 
method, we reduce consideration to only preimages of special points, i.e. points from 
a tiny cone- like set in the image. At his point our proof goes independently, because 
investigation of the number of regular preimages of special points is more concrete 
problem. 

Let us present a more detailed description of each step of the proof. 
1.2. Normal forms of local families and their integration. This step is done 



in [IY3] §0.3 and §1. We just say several words about it. 

It turns out that in a small neighborhood of an elementary equilibrium point 
there exists a finitely differentiablc normal coordinates (in the Cartesian product of 
the phase space and the parameter space) , so-called normal forms of an equilibrium 



point. The list of finitely differentiable normal forms was obtained in [ [Yl |. The 
main feature of the list: all normal forms are polynomial and integrable. The smaller 
the neighborhood of a normal form, the higher its smoothness. So smoothness can 
be chosen arbitrary large. All normal forms are summarized in Table 1 sect^. 

In a small neighborhood of an elementary equilibrium point one can choose two 
small segments, say £~ and E + , transversal to the vector field for the critical value 
of parameter and explicitly calculate the Poincare (correspondence) map which 
maps a point from one segment say £~ along the corresponding phase curve to a 
point from the other segment S + (see Fig.l). For an appropriate choice of segments 
£ — , E + and coordinate functions x, in E _ , S + respectively, and a smooth function 
A(e) in the original parameter e of the family the Poincare return map A e : x — > y 
can be explicitly computed. Moreover, there is a Pfaffian (with polynomial coeffi- 
cients) 1-form lo of the form 

(5) P{x, y, A(e)) dx + Q(x, y, A(e)) dy = 

which vanishes on the graph y = A e (x). For example, in the case of a nonresonant 
saddle A c (x) = x x ^ and lo = x dy + X(e)y dx. See Table 1 for the other cases. 

1.3. Singular-regular systems determining the number of limit cycles. 

We present a description of a system of equations determining the number of limit 



cycles. For a detailed description we refer to [ fY3 §0.4 and §1.4. 

Let 7 be a polycycle, occurring in a generic k parameter family, with equilibrium 
points pi, . . . ,p n (possibly with repetitions) and connecting phase curves 71, . . . , 7„ 
such that 7j connects equilibria pj with Pj+i respectively. For each 1 < j ' < n 
endow the point pj with a C r -normal coordinate charts Uj. Consider transversal 
segments "entrance" and "exit" which are parallel to coordinate axis of 
the normal chart. The phase curve 7j-i enters the neighborhood Uj through S~ 
and the phase curve jj exists Uj through £+. The normal coordinates induce 
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coordinates Xj and yj on Sj and E" 1- respectively. For some parameter values the 
corresponding vector field defines the following collection of Poincare maps: 

Aj(;e) : x 3 V] = A^x^e), j = l,...,n 
e ) : 2/j ->• = fi(yh e )i J = 1, • • • ,n ( mod n), 
where Aj(-, e) is a local Poincare map form the "entrance" segment E~ to the "exit" 
segment E+ and fj(-, e) is a semilocal Poincare map along the phase curve "fj form 
the "exit" segment E+ to the "entrance" segment EJ +1 . 

Now we decompose the monodromy map (the Poincare first return map) along 
the polycycle 7 into the chain of the local singular maps Aj and the semilocal 
regular maps fj of the total length 2n. Limit cycles correspond to the fixed points 
of the monodromy. But instead of writing one equation for the fixed points of 
the monodromy we consider a system of 2n equations, which will be called the 
preliminary basic system: 



(7) 



iy j =A j (x j ,e) j = l,...,n 
\xj+i = fj(yj,e), 3 = 1, • • • , n ( mod n) 

Recall that C r -normal coordinates on S ■ and yj 's are C r -normal coor- 

dinates on E^ . Thus the system involves C-smooth regular functions fj 's and the 
maps Aj from the list (modulo reparametrization e — > A(e)), that are essentially 
singular. The problem now is to estimate the number of solutions uniformly over 
all sufficiently small parameter values. 

1.4. The Khovansky reduction method. The system (]?]) is not easy to analyze, 
because it has the singular functions A,. The first key idea of the second step is 
to replace these singular equations in (0) by the Pfaffian (polynomial differential) 
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equations in the form (|5|). As a result we obtain the mixed functional-Pfaffian 
system of the form 

["j = 

^ \Fj(x,y,e) = j = l,...,n 

ujj = P-j dxj + Qj dyj, Fj(x, y, e) = x j+x - e) 

(x,y) = (x 1 ,y 1 ,...,x n ,y n ) G (R 2 ",0), ee (R fc ,0), 

where Wj are Pfaffian forms in the form (^) . This system can be interpreted as fol- 
lows: one has to take an integral manifold T for the Pfaffian equations of the system 
(||) and compute its intersection with the level set T~ x (0), where F : (R 2n , 0) — > R™ 
is the map with the coordinate functions Fj. In order to estimate number of iso- 
lated solutions to (]?]) one needs to estimate the number of isolated points in the 
intersection. It turns out that it is sufficient to analyze only transversal intersec- 
tions of r with a generic level set for b sufficiently close to the origin in 
R n . Since the integral manifold and the level sets have complimentary dimensions, 
a transversal intersection always consists of isolated points, which we call regular 
solutions to the system (@). What we are interested in is the upper estimate for 
their number, uniform over all the integral manifolds T and all sufficiently small 
values of the parameters. 



The method suggested by A. Khovanski [Kb] allows us to replace a mixed 
functional-Pfaffian system of the form (|J) by the two systems of a similar form, 
but containing n — 1 Pfaffian equations, n "simple" functional equations, and one 
special functional equation: the number of regular solutions to the initial equation 
is bounded from above by the sum of the number of regular solutions to these two 
auxiliary systems. 

1.5. ap-stratification and Bezout's theorem for a chain map P o F with 
a generic F. In this section we shall discuss the formula (@). The problem of 
estimating the maximal number of small isolated preimages is equally difficult for 
a chain map P o F : R n ->• R" with a generic map F : 1" -> R N , N > n and for 
a chain map P o j n F : R n — * R 71 with the n-jet of a generic map. We shall show 
that if the map F (resp. j n F) satisfies a transversality condition in an appropriate 
space, then F (resp. j n F) can be replaced by its linear part and we can apply 
the Bezout theorem to estimate the maximal number of small inverse images of 
the chain P o F (resp. P o j n F) uniformly over all sequences of numbers e\,...,e n 
decreasing sufficiently fast to 0. So, to simplify notations we shall consider a chain 
map of the form P o F : R n -> R". 

1.5.1. A Heuristic description. Consider a chain map P o F : R 2 — > R 2 , where 
F : R' 2 -> R N is a generic C k smooth map, k > 2 and P = (Pi, P 2 ) : R^ -> R 2 
is a polynomial of degree d. Fix a small positive r. We would like to estimate the 
maximal number of small preimages 



(9) ff{x G B P (0) : Pi o F(x) = e, P 2 o F(x) = 0} 

for a small enough e. 

To show the idea put N = 3, P\(x, y, z) = x 2 + y 2 , and P2(x, y,z) = xy. Assume 
also that F(0) = 0. Denote the level set by V e = {Pi = e, P 2 = 0}. The level set 
V e for e > consists of 4 parallel lines (see Figure 2). 
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Notice that in our notation the number of intersections of F(B r (0)) with V e 
equals the number of preimages of the point (e, 0) (^). 

It is easy to see from Figure 2 that if F is transversal to Vq it is transversal to 
V e for any small e > 0. Moreover, the number of intersections F(B r (0)) with V e 
equals 4 (see the points Pi, . . . , P4 in Figure 2). 

Another way to calculate the same number is as follows. Let us replace F by 
its linear part L F at zero. Then #{2; e B r (0) : Pi o F(x) = e, P 2 ° F(x) — 0} = 
#{x € B r (Q) : Pi oLp(x) — e, P2oLp{x) = 0} and solving this polynomial system 
also yields 4. 




Figure 3. The Idealistic Example 



The idea behind this picture is the following: Consider an arbitrary TV and a 
polynomial P = (Pi,P 2 ) : R N — > M 2 of degree at most d, N > 2. Define the 
semialgebraic variety V e = (Pi, P2)~ 1 (e, 0) as the level set. 

Assume for simplicity that for any small e / the level set V e is a manifold of 
codimension 2. We shall get rid of this assumption later (see Theorem ^b)). It 
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turns out that there exists a stratification of Vq by semialgebraic strata (Vo, Vq) (a 
decomposition of Vq into a disjoint union of semialgebraic sets see definition |30| ), 
depending on P only, such that 



(10) F is transversal to (Vq, Vq) 



F is transversal to V e 



Condition (Jlfj) is written for n = 2. Below we shall use its analogue for an arbitrary 
n. Let us present the key Proposition below and the simple of proof of it. This 
proof gives an insight to the main idea of the third step. 

Proposition 1. Let B r (a) be the r-ball centered at the point a£l 2 and let Lp, a 
denote the linearization of F at the point a. Under condition (f?7|), the number 
of intersections of the image F(B r (a)) with V e coincides with the number of in- 
tersections of the image Lp^ a (B r (a)) with V e , provided r is small enough. That 
is 

#{xGB r (Q): (P u P 2 )oF(x) = (e,0)} = 
( } #{i£B r (0): (Pi,P 2 )oL F)O (x) = (e,0)}. 

The argument below is independent of the codimension of V e . We only need 
condition ([l0]) and the fact that the codimension of V e coincides with the dimension 
of the preimage of a chain map P o F. 

Proof Consider the 1-parameter family of maps F t = tF + (1 — t)Lp deforming 
the linear part of F into F. Clearly, Fi = F and Fq = Lp. Fix a small r > 0. 



Since, F is transversal to Vq at all F t are transversal to Vq at 0. Condition (10) 
implies that for all small e and all t € [0, 1] F t is transversal to V e . 

Therefore, the number of intersections of Ft(B r (0)) with V e is independent oft. 
Indeed, assume that #{F tl (B r (0)) n V € } ^ #{F t . 2 (B r (0)) n V t } for some t x < t 2 . 
Then as t\ increases to t 2 there is a point t* where the number of intersections 
drops or jumps. At this point t* the condition of transversality of Ft* and V e must 
fail. This completes the proof of the proposition. 

2. Normal forms for local families and their applications. 

In this section we present the functional-Pfaffian system whose number of solu- 
tions bounds from above the number of limit cycles. This system was obtained in 
|IY§. 



2.1. Local families and polynomial normal forms. A local family of planar 
vector fields is the germ of a map, 

v : (K 2 , 0) x (R k , 0) -» (R 2 , 0), (x, y, e) ^ v(x, y, e). 

A C"*-smooth conjugacy between two local families v and w of the above form is a 
map 

H : (K 2 , 0) x (K fe , 0) -> (R 2 , 0), (x, y, e) » H(x, y, e), 

such that 

H*v(x, y, e) = w(H(x, y, e), e), 

where stands for the Jacobian matrix with respect to the variables x,y. (this 
definition does not yet allow for reparameterization of a local family) . Two families 
are finitely diffcrcntiably equivalent, if for any r < oo there exists a C r -conjugacy 
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between them. The two families v, w are orbitally equivalent, if there exists the 
germ of a nonvanishing function <fi : (R 2 , 0) x (R fc , 0) — > R 1 such that v is equivalent 
to (f) ■ w. 

To allow for a reparameterization of local families, we say that a family v(-, e) is 
induced from another family w(-, A), A G (R m , 0), if «(-, e) = w(-, A(e)), where A(e) 
is the germ of a smooth map (R fc , 0) — > (R m , 0). The number of new parameters m 
may be different from k. 

Assume that the family w(-,X) is global (i.e. the expression w(x,y, A) makes 
sense for all (x,y,\) G R m+2 ); this happens in particular when u; is polynomial 
in all its arguments. Restricting the parameters A onto a small neighborhood of a 
certain point (0,0, c) G R 2 x R m , we obtain a localization of the global family w, 
which formally becomes a local family after the parallel translation A ^ A - c. 



Definition 6. 1. A local family v = v(-, A) is finitely smooth orbital versal unfold- 
ing ( in short, versal unfolding) of the germ v(-, 0), if any other local family unfolding 
this germ is finitely differentiable orbitally equivalent to a family induced from v. 

2. A polynomial family w(-,X), A G R m , is a global finitely smooth orbital versal 
unfolding (in short, global versal unfolding) for a certain class of local families of 
vector fields, if any local family from this class is finitely differentiable orbitally 
equivalent to a local family induced from some localization ofw. 



To investigate a versal unfolding means to investigate at the same time all smooth 
local finite-parametric families which unfold the same germ v{-, 0). The main result 
describing versal unfoldings of germs of elementary singularities on the plane, is 
given by the following theorem. 



Theorem 7. |1Y3] Suppose that a generic finite-parameter family of smooth vector 
fields on the plane possesses an elementary singular point for a certain value of the 
parameters. If this point has at least one hyperbolic sector, than the family is finitely 
differentiable orbitally equivalent to a family induced from some localization of one 
of the families given in the second column of Table 1. 



Table 1. Unfolding of elementary equilibrium points on the plane. 
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Type 



Normal forms 



Poincare 
Correspondence maps 



Pfafnan equations 



So 



V = -Ay. 
A = A € K 1 



y = x, 
x > 0, y > 



x dy — Xy dx = 



x = x (i+P M ( u ,A)) 

y = -?/• 

u = u(x,y) =x m y n , 
P„{u,X) = ±m m (1 + A m m" 
+W„-i(it,A), 
A = (Ai, . . . , A M ) 



= mlogy + 

rv n du 

x > 0, y > 



y P M (y",A) - 
xP M (x m ,A) dy = 



D° 



p 



x = Qij.(x, A), 
y = -2/' 

Q M (x,A) =±x" +1 (l + A M x' 1 ) 
+W M _i(a;,A), 
A = (Ai, . . . , A M ) 



y = c(x)x 



x, y £ R 1 



x dy — y dx = 



= log 2/ + 
rl du 
■fx Q M (u,A) 

y > 0, x el 1 



M (x, A) dy 
y dx = 



In what follows the following notation for elementary equilibria (the subscript 
indicates the degree of degeneracy) : 

So — Nonresonant saddle; 

— Resonant saddle whose quotient equation (the differential equation for 
u = x m y n below) has the singular point of multiplicity fj,+ 1 at the origin, ji > 1; 
if we want to specify explicitly the resonance between the eigenvalues, we use the 
extended notation S^"™ 1 ^ assuming that the natural numbers m, n are mutually 
prime; 

Dp — Degenerate saddlenode of multiplicity n; 

W M _i(z, A) = A + X\z + • ■ ■ + A^-iz^ -1 is a Weierstrass polynomial of degree 
fi — 1. 



Different technical remarks concerning this table see in [[Y3 §1.1. We just briefly 
describe each column. 

The first two columns do not need extra words. In the third column of the 
table the Poincare correspondence maps y — A(x, A) for the polynomial normal 
forms are given. They are implicitly defined by the equations relating x to y, 
these equations depending explicitly on the parameters A and thus implicitly on 
the original parameters e. The choice of segments transversal to the phase curves 
of the family described in fig. 1. 



2.2. Basic system. Here we describe the system of equations which will be ana- 
lyzed from now on. Assume that a polycycle occurs in a generic fc-parameter family 
of vector fields, and all the vertices of the polycycle are elementary. 

Then the number n of vertices is < k. Moreover, one can claim that each vertex 
is of one of the types S^, fij > 0, or D^ c , /x, > 1, and — ^ ( see [IY2] §1.4). 

Next, we proceed with introducing the normalizing C p -smooth local coordinates 
near each elementary vertex, as this is described above (the exact order of smooth- 
ness will be specified later on). Then a pair of C p -smooth transversals may be 
chosen near each vertex, and endowed with local C p -smooth charts Xj , yj in such 
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a way that the correspondence map taking a point with a coordinate Xj on the 
"entrance" transversal to a point with the coordinate yj on the "exit" transversal, 
will be of one of the standard types listed in Table 1. 

More precisely, for each vertex j = 1 , . . . , n Theorem ^ yields the localization 
point Cj — (0, ...,0, Cj) G fJ 4 " 1 , where Cj G K 1 is the formal invariant of the 
unperturbed singular point, and also if j-th vertex is a resonant saddle, then the 
rational hyperbolicity ratio n : m is explicitly specified. 

Denote by Ai^(x, A) the correspondence map for each of the four types of singu- 
larities from Table 1, I = 5o, <Sp, or D^, with the corresponding index /i G N (for 
I = So by definition /i = 0). In case with /i > we consider the mutually prime 
pair of natural numbers as an additional parameter of the corresponding map, so 
in this case the rigorous notation would be As ^(x, A; [n,m]). 

Definition 8. 1. The unspecified basic system for determination of limit cycles 
occurring in k-parametric families of vector fields is the system of n regular and 
n singular functional equations in 2n variables Xj , yj , depending on parameters 
X j ,nj,mj,e, 

Vj = A l} iN { Xj , A^ ; [ nj , m 3 ] ) , A^ G + 1 , 
x j+x = fj (%,£), eG(M fe ,0). 

(12) j = 1, . . . ,n mod(n), h e {ScS^D^D*}, 

nj, rrij G N, G Z + , < k, n < k, 

A depends on rij,mj only if lj — with /i > 0. 

2. A specified basic system is one of a finite number of unspecified basic sys- 
tems together with an explicit indication of specification, which by definition is the 
collection of: 
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• localization points Cj — (0, . . . , 0, Cj) € { n particular this means that 
hyperbolicity ratios of all nonresonant saddles are explicitly given; 

• hyperbolicity ratios Uj : mj for all resonant saddles; 

• smooth functions fj(x,e) depending on the parameters e, are defined in some 
open neighborhoods (M. + ,0)j and fj(Q, 0) = 0; 

• characteristic size, that is, the value r > which determines the domain of 
the specified basic system as follows: 

0, y) G I r = {\xj\ < r, <r, j = l,...,n}c M 2 "; 
(A,e) eB r = {\\X-c j \\ <r, \\e\\ < r} C K *+/n+-+Mn , 

where X is the tuple of all parameters of all normal forms from Table 1, X — 
(A , ...,A n ); the characteristic size must be so small that all functions fj were 
defined for the corresponding values of their arguments. 

Notations related to definition |^ There is only a finite number of unspecified 
basic systems, each one being completely characterized by the string of discrete 
data 

(14) T = (h,fii, . . . ,Z„,Ai„) 

subject to the total restriction n < k, ^ fij < k. We call the data T the combina- 
torial type of the unspecified basic system. 
The string 

5 a = (c 1 ,...,c n ,...,m Ja ,ii 3t> ,...)er +2s 

(15) 

r > 0, Cj € K\ m 3a ,n 3a 6 N, fj G C p (M fc+1 , 0). 

will be referred to as the algebraic part of the specification (for reasons to be clarified 
later), while the string of functions 

Sf = f = (/!,...,/„) 

is called the functional part of the specification. The functions /j are defined on 
the domain I r x B r , where r is the characteristic size introduced earlier. 

Denote by B(T,S a ,f) the number of isolated solutions to the specified basic 
system (T, S a ,f) in the domain I ro . One can check that B(T,S a ,f) is defined in 
such a way that it bounds the cyclicity of the polycycle with such a specification. 

After all these notions (or rather the language) being introduced, we may for- 
mulate the problem of estimating cyclicity of elementary polycycles occurring in 
generic fc-parametric families as follows. 

Theorem 9. For any type T of unspecified basic system and any choice of the 
algebraic part S a one may choose the order of smoothness po and an open dense 
subset F = ^t, s a .r in the space of C Po -smooth functions C Po (I ro x B ro ,M. n ) such 
that for every f = . . . ,/n) G F and a sufficiently small characteristic size 
r'o = r (f) the number of isolated solutions B(T,S a ,f;r ) to the specified basic 
system (T, <S a ,f) in the domain I ro is uniformly bounded over all parameter values 
(A,e) G B ro : 



(lQp( T i S aJ;r()) = sup #{(x,y) satisfying (12), (x,y) £ I ro } < 2 

V ^ 0,A)£B ro 



25*- 



and, therefore, E{k) < 2 



25A-- 
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3. The Khovanski reduction method. 

In this section we describe the method of reducing a functional-Pfafhan system 
to a chain map of the form ([}]). The construction in its full generality is described 



in the book |Kh[ |. Our exposition relies on the one in |IY3], but has new important 



features so we can't just refer to neither |Kh|, nor [IY3] 



3.1. Pfaffian systems and their separating solutions. Let M be a smooth 
orientable n-dimensional manifold, not necessarily compact or connected, and to be 
a smooth 1-form on it. 

Definition 10. A codimension 1 smooth submanifold T C M is the separating 
solution for the Pfaffian equation uj = 0, if: 

a) T is the integral manifold, that is, the restriction of to on the tangent bundle 
of r is identically zero: 

Vx e r, Vv e T X T ui(v) = 0; 

b) r does not pass through singular points of lo : 

Vier, 3veT x M lv{v)\ TiM ^q- 

c) r is the boundary of a domain D C M and the coorientation induced on T by 
to, coincides with its coorientation as the boundary. In other words, on any vector 
pointing outward from D , the form is positive. 

Let now u)i, . . . , u>k be an ordered fc-tuple of smooth 1-forms on M. Consider the 
system of Pfaffian equations 

(17) wi = 0, ... ,w* = 0. 

Definition 11. A submanifold T is the separating solution for the system of Pfaf- 
fian equations, if there exists an increasing chain of smooth submanifolds, 



(18) L = T k C L fc _x . ••• 1; . 1:: M 

such that for any j — 1, . . . , k submanifold Tj is the separating solution for the 
Pfaffian equation on Tj^i, determined by the restriction of the form u>j on the 
latter submanifold. 

Let J- : M — > R s be a smooth map s < n — k. Recall that a point a £ M s is 
called a regular value for the map T if the linearization matrix, denoted by Jj?{x), 
has full rank for any x G J 7 ^ 1 (y). By the rank theorem the level set V a = J-^ 1 (a) 
of a regular value is a smooth manifold of dimension n — s. 

We call a e K s a regular value for T with respect to Pfaffian equations fiF\ ) if a 
is a regular value of F and the fc-form ft = lo\ A • • • A u>k, restricted to V a f2|y a is 
nondegenerate, i.e., singular points of tt\v a have measure zero. 

Consider a pair of smooth maps T : M -> W and F : M -»• K n - k ~ s . Now 



we add to a Pfaffian system (17) two types of functional equations. The first type 
consists of functional equations T = a, where o £ I s is a fixed regular value of 
T with respect to a Pfaffian system (p"^). The second type consists of functional 
equations F = b, where b 6 R ,l_fc_s is a variable. We call equations T = a, with 
a fixed a 6 K s , by rigid equations and F — b, with a varying b G IR ra_s_1 , by loose 
equations. 
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Definition 12. Let fl = (ui, . . . ,uik) 6 (A 1 (M)) k be a k-tuple of smooth 1-forms, 
T : M — > M. s and F : M —> R n ~ k ~ s be smooth maps, and a G R s be a regular 
value for T with respect to the k-tuple of smooth 1-forms. A solution to the mixed 
functional- Pfaffian system 

(19) fi = 0, F = b, T = a, beW l - k ~ s 

is a pair (T a ,Lb), where Lt, C M is the preimage F~ 1 (b) and T a is a separating 
solution for the Pfaffian system Q = 0, restricted to V a , and the intersection T a (~]Lb 
is nonempty. 

The solution is regular, if T a is the separating solution for the restriction of 
Pfaffian equations to V a and b is the regular value for the restriction of the map G 
on r a . If (r a , Lb)) is a regular solution, then the intersection T a n Lb is transversal 
and consists of isolated points. 

Definition 13. The Khovanski number JC{fl, F;J- = a} for the mixed system ( |7^ j 
is the upper bound for the cardinalities #{T a n Lb} over all regular solutions of the 
system. 

Remarks 1. The Khovanski number is also defined if k = (resp. s — 0), 
i.e., there are no Pfaffian (resp. rigid) equations at all. In this case one may put 
formally T = M (resp. V a = M), and JC{lo,F;J- = a} (resp. /C{f2, -F; 0},) is equal 
to the upper bound of the cardinality of preimages #{Lb H V a } (resp. #{Lb (~l T}) 
of regular values for the map G\v a ■ V a — > M. n ~ k ~ s . 

2. If we want to stress in the notation the phase space M of the functional- 
Pfaffian system, we use the notation K,m{Q,F,J- = a}. Usually this is necessary 
when F, T , and Q are defined on the Euclidean space W 1 , while we are interested 
only in solutions belonging to some (open) ball. 

3. If we fix a coordinate system in M. n ~ k ~ s , denote by F±,...,F S coordinate 
functions of the map F : M — > W l ~ k ~ s , and introduce the (n — k — s)-tuple of 
1- forms Qp = {dF\, . . . , dF s ), then we can consider the following mixed system 

(20) n = 0, Q F = 0, T = a. 

Regularity in the definition of the Khovanski number IC{Q,F; T = a} implies that 
lC{n, F; T = a} = IC{{Q, Q F ), 0; T = a}. 

The goal is using the Khovanski reduction principle estimate the Khovanski 
number for the mixed functional-Pfaffian system by a linear combination of the 
Khovanski number for some number of entirely rigid functional systems. 

The first step of the reduction principle is to estimate the Khovanski number 
for a given mixed system by a linear combination of the Khovanski numbers of two 
auxiliary systems containing a reduced by one number of Pfaffian equations and an 
increased by one number of rigid equations. 

The second step is using remark 3 replace all loose functional equations for 
pfaffian equations and apply the reduction principle to the mixed system consisting 
of (n— s) Pfaffian equations (f2, Qf) and s rigid equations. Thus, after (n — s) steps 
of the reduction principle we obtain a finite collection of entirely rigid functional 
systems. 
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3.2. The Reduction principle for one Pfafnan equation. We show how to 
eliminate the Pfaffian equation from the mixed system with [n — s — 1) loose equa- 
tions and s rigid functional equations. 

(21) lu = 0, F = b, T=a, F : M ~> R"~ s_1 , F:M-^W, 
We shall outline only the key ideas. 

Definition 14. A smooth positive function p : M — > R+ is called covering, if it 
tends to zero along any nonaccumulating sequence of points in M . In other terms, p 
vanishes "at infinity" on M, so that all level hypersurfaces of the covering function 
are compact subsets of M. 



Remark 1 . This definition applies both to compact and noncompact manifolds, but 
in the compact case a smooth function is covering if and only if it is everywhere 
positive, thus automatically bounded away from zero. 

Suppose that the manifold M is endowed with the Riemann volume. Since it is 
orientable, one may use the duality between functions and n-forms on M . Denote 
by the asterisk the operator taking an n-form into the function (dividing by the 
volume form). 

Fix Euclidean structures in R n_s_1 and R s . Let Fx, . . . , F„_ s _i and T\, . . . , T s 
be the coordinate functions of the maps F and T in (Klh respectively. 



Definition 15. The contact function for the mixed system (|£ip is 

(22) T s+1 = *(uj A dF\ A • ■ ■ A dF n _ s _ x A dT\ A • ■ ■ A dT s ). 

The operator taking the mixed system (lu, F; F) into the corresponding contact func- 
tion, will be denoted by a : (w,F\T) i— > a(uj, F; F) = J- s +i- 

Define the two maps by their coordinate functions, 

(23) T c = {Ti, . . . , F Sl Ts+i), J 700 = (Fi, ■ ■ ■ , F s ,p), 

both taking M to R s+1 , where ,F s +i is the contact function (|22|), and p is the 
covering function. 

Theorem 16. Suppose that the system admits regular solutions in the sense 
of Definition [7j]. Then for any sufficiently small regular e 

(24) K{w, F;J r = a}< hc{w, F; F°° = (a, e)} + JC{w, F; T c = (a, e)}, 

where regularity of e means that (a, e) is a regular value for both J-°° and T c and 
is necessary to the right-hand side systems being well defined. 

Before proving this theorem recall the Rolle lemma from an elementary calculus. 

Lemma 1. Consider C 2 Morse functions f : S 1 — > M 1 on the circle and g : [0, 1] — > 
R 1 on the segment, i.e., functions f and g have only nondegenerate critical points. 
Then for all a £K 

f(x) = a} < #{x : f'(x) = e} 
U '" 9{x) - a} < #{x : g'(x)=e} + l 

for any sufficiently small e. 
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Proof Prove the formula for / : S 1 — > M 1 . For a sufficiently small e the number of 
local maxima and minima equals #{x : f'(x) = e}. Between any two consecutive 
preimages x\ and X2 of a point a, i.e., f(x\) = /(a^) = a there exists a local 
minimum or maximum. Q.E.D. 

Formula ( |25| ) in the case of one equation transfers a loose equation into a rigid 
one. 

Proof of theorem Take a reg ular solution (r Q , L b ) for (gf|), where L b = F -1 (6), 
and suppose that the intersection T a n L& consists of isolated, say d, points. Since, 
b is regular value of the restriction F\ra , any small variation of b may only increase 
the number of intersections. Take b to be a regular value of the restriction F\y a 
or equivalently (b, a) to be a regular value of the map (F, J-) (rather than of the 
restriction of T to r a ). 

Then any level set Lb is a one dimensional smooth manifold, intersecting T 
transversally. By the classification theorem for one-dimensional manifolds, Lb is 
the union of compact (circles) and noncompact (lines) components. Fix some ori- 
entation on each circle and each curve in Lb. Consider the function /' : Lb — > K 
which maps a point x E L a to the value of the f-form lo on the unit positively 
oriented vector tangent to Lb at point x. 

Fix a connected component, denoted by 7 C !(,. Between any two consecutive 
intersection x and y of L a with T values f'(x) and f'(y) must have different signs. 
Now we can apply the Rolle lemma with /' = /', when 7 is a circle, and /' = g' , 
when 7 is a line. 

Each point x where f'(x) = (resp. f'(x) is small) is the point where the 
linear functionals dFi(x), . . . , dF n _ s _fe(x), dJ-'i(x), . . . , dJ- s (x). and u>(x) are linear 
dependent (resp. almost dependent), i.e. J- s +i(x) = (resp. !F a +i(x) — e). This 
completes the proof of the theorem. Q.E.D. 

Corollary 3. If the manifold M is compact, then for any sufficiently small regular 
e 

(26) K.{u,F;J r = a}<1C{u,F;F c =(a,e)}, 
where regularity of e means that (a, e) is a regular value for T c . 

Proof Indeed, in this case the first term in ( pl| ) disappears. 

Remark 2. The choice of the Riemann volume form is not essential for the above 
construction. Indeed, if the volume form vol n is replaced by a new one b-vol n , where 
b is a positive function, then the function !F s +x will be replaced by !F s +i = b^ 1 J- s +i, 
and the map T c = . . . ,T a ,T s +\) will have the same zero set. 

3.3. The Khovanski reduction in the general case. Consider now the general 
case of the mixed system ( p^ ) with k > 1. Suppose that T a is a separating solution 
for the Pfaffian system Q = restricted to V a . By definition, this means that 
there exists a separating solution Vj, C V a to the Pfaffian equation uik = on a 
separating solution r^_ x C V a to the Pfaffian system fl' = restricted to V a , where 
f2' = (ct>i, . . . , Note that if p is a covering function on the manifold M, then 

its restriction on T c j._ l is the covering function for the latter submanifold. Next, one 
can endow V a (resp. r^_ x ) by the Riemann (n — s)-volume (resp. (n — k — s + 1)- 
volume) form voly~ s (resp. vol^ k ~ s+1 ) in such a way that 

a k — 1 

(27) dT\ A • • • A dT s A vol^T* = vol^ ui A • • • A u> k -i A volp^ 1 = vol^T* 



18 



Bifurcation of polycyclcs 



Since the forms ujj, j = 1, . . . , k — 1 are linear independent in a neighborhood of 
rj-ij these formulas define volume forms near V a and respectively. As this 

was mentioned before, the choice of the Riemann volume form does not affect the 
assertion of Theorem |l7|. 

Thus one can apply Theorem [l^ to the mixed system 

(28) u k = 0, F=b, T = a 

on the manifold r°_ fe C V a . To describe the result, we introduce the following two 
maps from M to M"- fe+1 , 

(29) T C = (T U ...,T S , P ), = (F 1 ,..., T a ,T„), 

where p is the covering function on the manifold M, and : M ~ * R is the smooth 
function obtained as 

= <j{tt, F; T) = *(dF 1 A ... A dF n _ k _ s A 
AdTx A . . . AdT s A uj 1 A . . . A u k ). 

The above choice of the Riemann volume on r^_ 1 implies that the asterisk 
operator in the ambient manifold M agrees with the asterisk operator relevant to 



r|_ 1 , therefore the formula (30) defines the same function as the formula fl22|): 
LOk A dFi A • • • A rfF„-fe- s A dT\ A ■ ■ ■ A dT s = ■ volp7 fc+1 . 

' 1 k - 1 

Theorem 17. Let £l,F, T c , and J-°° be as above. Then for any sufficiently small 
regular e 

(31) £{fi, F; T = a} < \K{Q.' , F; T°° = (a, e)} + JC{VL\ F; T c = (a, e)}, 
where regularity of e means that (a, e) is a regular value for both F 00 and T c . 

Corollary 4. If either V a is compact or the restriction F\y a : V a — > M. n ~ k ~ s is a 
proper map, i. e. preimage of any point is compact, then for any sufficiently small 
regular e 

(32) £{Q, F; T = a} < K{tf, F; T c = (a, e)}, 
where regularity of e means that (a, e) is a regular value for T c . 

Proof Straightforward application of Theorem 

Rcrating the above two statements, one can replace one by one the Pfafiian 
equations by the rigid functional ones, obtaining new systems whose Khovanski 
numbers estimate from above that of the initial one, by virtue of the inequalities 
( pl| ) and its compact counterpart (|24|) . On each step one has two possibilities, either 
to replace a Pfafiian equation by the contact function, or by the covering function. 
But once the covering function appears among the rigid functional equations, the 
level sets F _1 (-) n V a becomes compact as a submanifold of a compact V a , hence on 
the next steps the Corollary to Theorem [IT] applies rather than the Theorem itself. 

Denote by T e c and T e °° the two operators, transforming the mixed system 
{n, F;T = a} into the mixed systems {fi', F; T c = (a, e)} and {Q', F, J 700 = (a, e)} 
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respectively, where the maps T c and T°° are given by ( p9| ) and (p0|): 

T e c {fi, F;T = a} = {tf, F; (T, F, F}) = (a, e))}, 

T~{n, F-T = a)} = {Of, F; (T, p) = (a, e))}. 

If we start with the mixed functional-Pfaffian system {(f2, SIf), w\J- = a}, with 
(f2, Qp) being an (n — s)-tuple (u>i, . . . ,u>k, dF±, . . . , dF n -k-s), and eliminate sub- 
sequently the forms u>k, w*— i, ■ ■ • , Wi, d-Pij ■ ■ ■ , dF n ^k-s, then the following maps 
from Af to R™ arise: 

a) the map J 7 ™ , if on each step the contact function was used, 

(34) { Wj ^ [0] - (a, e"- s ))} - (T e c _ o • ■ • c T e c J{0, F\T — a}, 

where e"~ s = (e 1( .. . , e n _ s ); 

b) the maps Tu\ , if on the jth step the covering function was used, while on all 
other steps the contact ones were, j = 1, . . . , n — s, 

m to = (*, ^ s ))} = PZ_ t o • • • o t* +1 o r- o 

oT £ c j _ i o...oT e c i ){f7,F;^ = a }. 

Then inductive application of Theorem [l6] immediately yields the following fun- 
damental result. 



Theorem 18. The Khovanski number for the mixed system on a manifold 
M with the covering function p and any sufficiently fast decaying to zero sequence 
(ei, . . . , e n -s) admits the upper estimate by a linear combination of Khovanski num- 
bers of some (n — s + 1) auxiliary systems, each of them containing only rigid 
equations and no Pfaffian equations at all: 

k 

(36) JC{Sl, F-F = a}< JC{lu; T n = (a, e"" 3 )} + \ £ >C{lu; = (a, e"" 5 )}, 



where the maps Tu\ are defined by the formula ft33[)-(p- 

3.4. Applications. The Khovanski reduction process is constructive. This leads 
to the result, which will be now formulated. 

Assume that the manifold M is an open domain in R™ and admits a polynomial 
covering function p. The main example is the unit ball B = {x 6 R" : x? < l}, 
for which one may take p(x) = 1 — ■ x?. Then the Riemann volume form can be 
chosen algebraic, dx\ A • • • A dx n . 

Assume also that all the forms u>i, i = 1, . . . , k are polynomial (i.e. with poly- 
nomial coefficients), and the maps T and G are at least C ra_s -smooth. Then, since 
the operators Tjf and T e °° introduced above, involve only algebraic operations and 
differentiation of functions, the following holds. 



Theorem 19. If the system has no rigid functional equations at all (s = 

0) and is defined on a semialgebraic subset M C R", all Pfaffian forms and the 
covering function p are polynomial of degrees < d, then all the maps J-\ a ] : M — > R n , 
a = 0, 1, . . . , k constructed in Theorem [7^ are of the form 

(37) T [a] =P a oj n - s F, 
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where j n s F is the (n — s)-jet extension of F, and P a are certain polynomials 
defined on the jet space J"- s (Jfr,R n - fc ). For all a = 0, 1, . .., k the degrees of the 
polynomial P a admits the upper estimate by 2 a (dk + n) and each map has a 
regular point for a generic map F . 

Proof The reduction procedure of elimination of a Pfaffian equation boils down 
to consecutive application (n — s) times of one of the operators T e c or T £ °° (Q). 
Consider the first step. 

Tl' a = a{(Q, Q F ), 0; F)} = A • • • A u k A dF 1 A • • • A dF n _ k ) = P 1 ^ o j 1 F, 

where P 1,a : J 1 (W 1 ,M. n ~ k ) — > R is a polynomial of degree at most dk + n and is 
defined on the space of 1-jets J^R™, R n - fc ). 

Denote by f2* = (0,0^) the (n — s)-tuple of the 1-form, dF s by lo^+s for s = 
1, . . . ,n — s, and the (n — s — r)-tuple of the 1-form, which consists of all of 1- 
forms of f2* except of the first r, by Q*. Consider T r * a = a{£l*_ 1 ,uj r ; J- r ~ 1 ' a } and 
T r ' a = (J 7r ~ i ' a , J^' a ) for r = 1, . . . ,n — s. It is easy to see that Tl' 01 has the form 
Fl' a = p r > a o j r F. 

Using induction in r it is easy to see that for r ^ a the degrees of correspond- 
ing polynomials p r ^ 1 < a and P r ' a , defined above, satisfy the following inequality 
degP r ' Q < 2 deg P r ~ 1,Q . For r — a the operator T £ °° will not exceed the degree 
dk + n and T r - a = p. This implies that degP r,a < 2 a (dk + n) and complete the 
proof. 



4. Functional-Pfaffian system for limit cycles 

In this section we consider a specified basic system (T, 5 a ,f;r) obtained from 
the unspecified basic system (g), that is we consider a system (g) together with a 
collection of formal invariants (ci, . . . , c„) of all singularities (which determines a 
point in the A-space) , a collection of hyperbolicity ratios nj a : uij a of all resonant 
saddles and a tuple of sufficiently smooth functions fj , on a sufficiently small open 
cube I r x B r in the (e, A)-space. 

Our local goal is to reduce this system to a functional-Pfaffian system having 
the form described in section with the following properties: 

• the new system has the form allowing for application of Theorem [l?]; 

• the number of regular solutions to the functional-Pfaffian system is greater or 
equal to the number of isolated solutions to (||), up to k, where k is the number of 
parameters of the original family. 

After application of Theorem [l^ we will obtain a number of chain maps with 
controlled degrees of the exterior polynomial parts. 

4.1. Upper estimate of the number of solutions for the basic system: 
statement of results. First of all we make the following remark. The algebraic 
part of the specialization can be identified with a point 

(38) S a = (ci, . . . ,c n ,n h ,m jl i . . . ,n 0s ,m js ) e R" +2s , 

where s < n is the number of resonant saddles on the polycycle: the fact that the 
numbers rij a , rrij a are in fact natural will become inessential for our constructions. 

Theorem 20. (reduction from basic to functional-Pfaffian system) Consider an 
unspecified basic system of a certain type T in codimension k, together with an 
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arbitrary specification 

S=(S a ,f;r), 5„el" +2s , f e C p {I r x B r , R n ), r > 0. 
Then one can explicitly construct a functional- Pfaffian system of the form {Q, F}, 
CI = (Hi, . . . , Q n +2s), F = ■ ■ ■ , F n+ k+m), m — X) I 1 ]; or * n a m ore traditional 
notation, the mixed system of loose functional and Pfaffian equations (no rigid 
equations) 

(39) A = 0, F = a; 

defined in a certain open bounded semialgebraic subset 

M = M(r) C I r x B r x R 2s 

(see Definition |8|), such that the following holds: 

• For any choice of the parameters (e, A) G B r the number of isolated (x,y)- 
solutions, denoted by B(T,S a ,f;r) of the specified basic system (T, S a , f ; r) admits 
the estimate by the Khovanski system ( p^ ) on the manifold M = M(r): 

B{T,S a ,i-r) <lC M(r) {n,F;Lj} + k; 

• The forms ilk have coefficients which are polynomial in all their arguments, and 
also in coordinates of the point S a G M" +2;i ; the degrees of those polynomials do not 
exceed 6fi + I, where fi is the order of degeneracy of the corresponding equilibrium 
point; 

• The covering function p(- ; r) for the phase space M(r) is polynomial in all its 
arguments and also in r, of the total degree not exceeding lAk; 

• The coordinate functions of the maps Fp are explicitly given as polynomials of 
the first degree on the 0-jet space of functions J°{I r x £? r ,R") with coefficients ±1. 



The proof of this theorem is completely construc tive and given in [IYS]. We only 
point out degree estimates which are not given in I Y3 1 . 



22 



Bifurcation of polycycles 



Table 2. Separating solutions for Pfaffian systems associated with unfolding of 
elementary equilibrium points. 



Type 


Submanifold 7 


Domain M r , Covering function p 


Piarhan system il = 


So 


y = A(x,X) 


< x, y < r, 
X G L r , 
p = xy(r - x)(r - y)p 


x dy — Xy dx = 


S M 


y = A(x,A) 

2 =X m 
n 

w = y 


< x, y, z, w < r, 

X £ L r , 

P„(z,X)^0, 
p = xyzw(r — x)(r — y) x 
(r - z)(r - w)Pp(x, X)p 


x dz — m z dx = 0, (1) 
y dw — n w dy = (J, (2) 
mP M (w, A) x 
yPfj,(z, A) dx — 
(mP^(w, A) +n) x 
i^(2,A)<i!/ = (3) 


D% 


y = A(x,X) 


\x\, \y\ <r, x^O, 
p=^--x%{?-y^~p 


x (x dy — y dx) — 




y = A(x,X) 


<y < r, \x\ < r, 
Xe L r 

Qn(;*)\[x,i) > 0, 
p = y(r - y)(r 2 - x 2 )x 

Q^{x,X)p 


Q ll (x,X) dy — y dx —0 



Notes to the Table Here we use the same notation as in Table 1 (and in fact 
this Table continues Table 1). In particular, n : m is the hyperbolicity ratio in the 
resonant saddle case S^. 

In the third column of the Table the symbol L r stands for a small r-cube in the 
(/i + l)-dimensional space of the parameters A, centered at the localization point 
c = (0, . . . , 0, c) G R A ' +1 , corresponding to the unperturbed system: 

L r = {Aer +1 : I Aij <r, i = 0, . . . , fi - 1, |A M - c\ < r} 

Everywhere in the Table the function p = p(A) is the covering function for the set 
L r , defined as 

p^^^-xD-.-^-xl^-^-ic-x^). 

This is a polynomial of degree 2/z in all variables A, r, c. Recall that degP m u = 
2p, and degQ^ = 2p + 1 (see Table 1). Each covering function p is therefore a 
polynomial (explicitly written in the Table). Thus, we obtain the following degree 
estimates: 

Type So: degil = 1 and deg p = Ap + 4. 

Type S^, p > 0: degfl < 6p + 1 and degp = 6p + 8. 

Type D^: deg ft = 2 and degp = 2p + 6. 

Type D^: degfi = 2p + 1 and degp = 4p + 5. 

Along with the estimate ^ & (the sum of codimension) this gives the 

estimates deg 17 < 6p, + 1 and degp < 14fc. 

4.2. Principal functional— Pfaffian system. We proceed with writing down the 
principal functional-Pfafhan system explicitly. Slightly abusing notation, we add 
the subscript j for objects related to the jth singularity, while letters without this 
subscript refer to objects related to the entire polycycle. In this notation we omit 
the reference to the characteristic size, still keeping in mind that all formulas are 
explicitly polynomial in r. 
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Notations Denote by Mj the domain from Table 2, associated with the j-th 
singular standard map, let jj be the corresponding manifold (separating solution) 
and by Q.j the tuple of Pfaffian forms on it: if the singularity is of the type D M or 
So, then £lj consists of only one form ujj — Aj dxj + Bj dyj, while in the case S^, 
p > 0, there are three forms, of which we denote the third one by Aj dxj + Bj dyj, 
(see Table 2). The covering function for Mj is denoted by pj : Mj — ► R+. 

Construction of the principal system The phase space for the principal functional- 
Pfaffian system is the Cartesian product of phase spaces corresponding to all the 
vertices of the polycycle and the r-cube in the e-space: 

M = M (r) = Mi x ■ ■ • x M„ x B r , 

(40) Mj = Mj ir are taken from the second column of Table 2, 

B r = {M <r,i = l,...,k}. 

Dimension of the phase space is equal to 2n + 2s + k + m, where: 

• k is the number of the parameters e (the principal integer index) ; 

• n < k is the number of vertices; 

• s < n is the number of resonant saddles on the polycycle (each such a vertex 
contributes two additional variables Zj , Wj into the list of independent variables) ; 

• m — n + f^j < n + k < 2k is the number of additional free parameters 
A= (A 1 ,..., A"), A^ el"i +1 . 

The covering function for such a space is the product 

(41) p = P f-p n -p e :M^R 1 + , 

where the last factor is the covering function for B r . From Table 2 it is clear that 
p is a polynomial of degree at most J^ji^^j + 8) < 14fc in both phase variables and 
the characteristic size r. 

Each form on Mj can be pulled back on M , yielding the form which is indepen- 
dent of all the coordinates except for those related to the jth. vertex. Denote by Q, 
the union of the tuples flW; thus is itself the tuple of 1-forms on M, containing 
n + 2s of them: 

n = (n (1 \...,n^) = (n 1 ,...,n n+2s ), 




{ujj} if j is not a resonant saddle, 
{ojji, u)j2, Wj} otherwise, 



where u>j\ = rrijXj dzj — Zj dxj, Uj2 — n-jyj dwj — Wj dyj, 

u)j = Ajdxj+Bjdyj. 

Each 7j is a separating solution to the Pfaffian equation or system of equations 
= on Mj, therefore the Cartesian product 

T = 7i x ■ • • x 7„ x B r 

is the separating solution to the Pfaffian system f2 = on M. Indeed, one may 
consider the chain of submanifolds 

Ti = 7i x • • • x 7,; x M i+ i x • • • X M n X B r 

This chain possesses all the properties required by the definition of a separating 
solution, see section ||: there are no singular points of Pfaffian forms on all the 
manifolds from this chain, and the topological condition of Ti+i being the boundary 
of a domain in Ti is trivially satisfied, because each 7 i+i is the boundary of the 
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corresponding subdomain in Mj . Thus the Pfaffian part of the principal system is 
constructed. 

In this Pfaffian part we have the following information about the polynomials 
(recall that S a stands for the algebraic part of the specification for the basic system, 
which is identified by (|3|) with a tuple of real variables) : 

Aj,Bj 6 Z[x, y, A, S a ], degAj,degBj<Qp,+ l, 

p G Z[x, y, A, e,r], degp<14k. 

Now we proceed with description of the functional part of the principal system. 
It is given by the map 

F = (Fi, F n+k+m ) : M - K»+ fc +™ 



(44) 



fj(.Vj,t)i j = l,..., n mod (n), 
Fj = {ej- n , j = n + 1, . . . ,n+ k, 

Aj_ n _fe, j = n + k + l,...,n + k + m. 

The dimension of a generic fiber is equal to the codimension of separating 

solutions of the Pfaffian system. An essential feature of the above map is the 
following one: the coordinate functions of the map F are polynomial combinations 
of the coordinates on the source space and generic functions fj : 



(45) Fj G Z[x, e, A, f], degFj = 1, 

and all coefficients of those polynomials are ±1. A more invariant way of formulat- 
ing the same property is to say that F is a polynomial map defined on the space 
of 0-jets of vector-functions 

(46) f : M — ► R n , 

and this phrase makes sense since M is a subset of a Euclidean space. 



Definition 21. The functional-Pfaffan system with the Pfaffian equations 
the functional equations defined on the domain ( f^| ) considered with the cov- 

ering function (j^/j), will be called the principal functional-Pfaffian system. The 
information provided by the estimates allows us to say that the principal 

system is effectively described. 



Later on we will refer to the principal system as simply the system (p9|). 

4.3. Reduction to singularity theory. The system (|39|), whose Khovanski num- 
ber majorizes the number of solutions to the basic system ([l2|), satisfies the con- 
ditions of Theorem |l8|. The conclusion of the latter claims that the number 
/C{f2, F;uj} is in turn majorized by the combination of Khovanski numbers for 
some 2n + 2s + fc + ro+l entirely rigid systems (recall that n + 2s is the number 
of Pfaffian equations and n + k + m is the number of loose functional equations in 
the principal system, which should be eliminated). The properties of the principal 
system, listed in the formulation of Theorem ^0], yield a complete description of the 
resulting systems as chain maps (the definition is given below). 

In what follows we treat the original variables Xj , yj , the auxiliary variables 
Zj a , uij a and the parameters e, A in almost the similar way, as it is suggested by the 
functional equations (|||) of the principal system (^). The algebraic part S a of 
the specification, however, plays a different role: the coordinates of the localization 
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points cj and the integers rij a , m 3a determining the hyperbolicity ratios of resonant 
saddles, would determine the point in the new phase space, around which the 
resulting chain maps will be considered. Recall that in §1 we introduced the vectors 
cj and c as 

cj = (0, . . . , 0, cj) G W*i +1 , Cj e M. 1 , 
c = (ci,...,c, i )et m , m = n + y^/ij. 

For our purposes it would be convenient to consider all (new) variables as taking 
values around the origin in the corresponding phase space. For this sake we make 
a parallel translation in the A-space, which would take the origin into the point 
c. Clearly, this translation does not affect the algebraic structure of the principal 
system (p9|), though changes the appearance of the equations. 
The characteristic size r retains its original meaning. 

Notations According to what has been said, we introduce the following nota- 
tions: 

X = {x,y,Z,W,e,\-c) 6 R2»+2«+*+*n 

f = (fl, ■■■Jn), fj=fj{Vj,£) f = f (x), 

where f is now considered as a vector-function of the argument x, though each 
coordinate function fj of the vector f depends in fact only on some of the coordi- 
nates of the vector x. By Z? p f we denote the collection of all partial derivatives of 
functions fj of the order p. 

We will also use the same notation M(r) for the domain of the principal system, 
though in fact it would become a subset of the unit cube ||x|| < r centered at the 
origin in the x-space. 

Now we can formulate the properties of the systems of equations which appear 
after elimination of Pfafhan equations from the principal system (p9[) as this was 



described in S 3.4 . Let m = 2n + 2s + k 



Theorem 22. Let m = In + 2s + k + m. For any fixed combinatorial type T 
of the principal functional- Pfaffian system any choice of the algebraic part 

S a of the specification and sufficiently fast decaying to zero sequence of numbers 
ei,...,e m , the number of nondegenerate solutions to the principal system in the 
domain M(r) for any choice of the characteristic size r > does not exceed the 
sum of the Khovanski numbers for m + 1 entirely rigid system of equations in the 
same domain. Each of these systems has the form 

(47)P(x,f(x),^ 1 f(x),...,^ m f(x);5 Q ,r) = (ei, . . e m ), x e M(r) C R m , 
where 

• m < 7k is the total number of variables ( the dimension of the phase space); 

• P is a vector polynomial, P = (Pi, . . . , P m ), Pi £ Z[x, . . . ;<S a ,r]; the degrees 
of each polynomial Pi is bounded by lAk2 l i = 1, . . . , m; 

• the domain M(r) belongs to the r-cube of the space R m , centered at the origin. 

4.4. Chain maps and related finiteness theorems. Now we proceed with a 
more invariant description of the geometric object corresponding to the system of 
equations (Em 
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Definition 23. Let K m be a Euclidean space with a fixed coordinate system x = 
(Xi, . . . , X m ), and U C R m a domain of the rectangular form, 

U = {a, < Xt < (3 l , i = l,...,m}. 

Denote by I the index subset I = {1, . . . , m} enumerating the coordinates in R m , 
and let for any j = 1, . . . , n Ij be a nonempty subset of I, 

0jHjCI, j = l,...,n. 

We say that a vector-valued function 

f:[/^K n , f =(/i,...,/„), 

is a Cartesian function of the Cartesian type X — (7i, . . . , I n ), if for any j the jth 
component of this function depends only on the coordinates Xi with i e Ij : in other 
words, 

For any given Cartesian type X with n = 1 the set of all C p -smooth Cartesian 
functions (ieCartesian maps with n = 1) of this type constitutes a Banach space 
with the natural C p -norm. We denote this space by Cj, sometimes omitting the 
explicit reference to the type X when the latter is clear from context. The space Cj 
will be referred to as the Cartesian space. In the same way the Cartesian spaces 
of maps arise. As a consequence, we may say about genericity of Cartesian maps 
(functions) within the given Cartesian type; the notions of openness and density of 
subsets are also naturally defined. 

Definition 24. Let f be a C p -smooth Cartesian map of a given Cartesian type X, 
and s > an nonnegative integer number, s < p. A Cartesian s-jet of the function 
f at a point xo G U is the equivalence class of all Cartesian functions of the same 
Cartesian type, which differ from f by a term which is s-flat at xo-' 

j s f(x ) = {geC p :|f-g|=o(|x-xo| s )}. 

The space of all s-jets of functions of the given Cartesian type X at all points 
xo £ U will be denoted by Jj(M m ,]R n ) or simply by 3 s , when the environment is 
unambiguously defined by the context. 
The map 

x^ff(x) 

is called the Cartesian s-jet extension of the Cartesian map f . 

The space of Cartesian jets of any type and any finite order admits a natural 
coordinate system, in which the Cartesian jet extension of a map f = (/i, . . . , f n ) 
takes the form 

x =(x 1 ,...,x M )^ (x,.F(X),{!^, teij 

all partial derivatives of functions fj of all orders up to~) 

s in the variables on which each fj actually depends J 

The Cartesian jet spaces possess almost all properties of the standard jet spaces. 
In particular, the natural projections 

(48) RM D V , y„ j0 ^ RM x rK t pr, _ _ _ | pr. JS (P r. + 1 _ _ _ 
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are well defined and endow each JJ- with the structure of an afHne bundle over 
R m . Thus it makes sense to say about polynomial functions defined on Cartesian 
bundles. 

Definition 25. A chain map with the exterior part P and the interior part f is a 
map of the form 

R m D[/9xn P(j|f(x)) G R m , 

where: 

• f is a Cartesian map from a certain Cartesian space Cj(R m ,R n ), and jj is 
the corresponding s-jet extension of f ; 

• P : JJ(R m ,R n ) — > R m is a vector polynomial (eventually depending polynomi- 
ally on some additional parameters), 

• the composite map is between the spaces of the same dimension:dirax = 
dim P = m. 

Having introduced the notions of Cartesian functions, maps, jets etc, we can 
describe the system (3.10) as a chain map defined on a small cube of some size 
r > with the exterior part P which is a polynomial with integer coefficients 
and of a controlled complexity; this polynomial depends on r and some additional 
variables A as well, and the interior part f belongs to some Cartesian space, since 
the functions fj depend only on some components of the vector x = (x, y, z, w, e, A — 
c) (recall that all nonzero coordinates of the vector c are already included among 
the variables S a ). Thus our problem of estimating cyclicity of a polycycle takes the 
following form: describe the Cartesian maps f for which the chain map admits an 
upper estimate for the number of preimages of regular values. 

Consider chain maps of the form 

x^G r (x)=P(j|f(x),r) = (P 1 ,...,P m )(j|f(x),r), x e U C R m , r > 0, 
depending polynomially on an additional variable r, so that 

(49) P : Jf(R m ,R n ) x M 1 -> R m , feC^f/,1"). 

We assume that the polynomial P and the Cartesian type I are fixed (and U 
denotes as before a unit cube) and P is nontrivial polynomial, i.e. at some point 
x G U the linearization matrix dP{x) has full rank. 

Suppose that the smoothness order p is sufficiently high, 

p > m + 1. 

Theorem 26. For any polynomial P = (Pi, . . . , P m ) as in fcfy) one may choose a 
subset Fp C Cj(U, R n ) in the space of Cartesian functions of the given type, which 
is open and dense in this space such that for any Cartesian function f £ Fp and 
any sufficiently quickly decaying sequence a%, . . . ,a m there exists a characteristic 
size ro > such that the number of preimages of {a\, . . . , a m ) admits the following 
upper estimate: 

m 

(50) limsup#{x x 6 [/ n G r (x) = ( ai ,...,a m )} < TTdegP, < 2 25fc2 . 

A bit of terminology: "Replace an n-th jet j n F by its linear part at a point 
a e R n " means "replace the map j n F : R™ -> J"(R",R n ) by its linear part L F , a ,n 
at the point a" . 
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By the phrase "a map G : M — > N of manifolds satisfies a transversality condi- 
tion" we mean that for some manifold (resp. a collection of manifolds) in the image 
N the map G is transversal to this manifold (resp. these manifolds). 

The second stage consists in constructing a stratification of the n-jet space 
J n (R",R") (a decomposition into a disjoint union of manifolds described below) 
such that if the n-jet j n F is transversal to all manifolds of this stratification, then 
the following theorem is true: 

Theorem 27. Let P = (P\ , . . . , P n ) be a nontrivial polynomial defined on the space 
ofn-jets P : J n (R n ,R n ) -> R" and let F : R" -> R" be a C k smooth map, k > n. 
Suppose the n-jet j n F satisfies a transversality condition depending only on P. 
Then for a sufficiently small r one can replace in the statement of the previous 
theorem the n-jet j n F at the point a by its linear part Lp,a,n- Namely, 

#{xeB r (a): P 1 oj n F(x) = a 1 ,...,P n oj n F(x)=a n } = 

#{x e B r (a) : PioL Fta>n (x)=ai,...,P n oL F>a>n (x) = a n }, 

where a\ , . . . , a n go to zero sufficiently fast. By Bezout 's theorem the number of 
solutions to the equation in the right-hand side of can be bounded by the product 

rnudegp,. 



The classical transversality theorem [AGV] says that for a generic map F its 



n-jet j n F satisfies any ahead given transversality condition. 

4.5. Stratified manifolds. Now we recall basic definitions from the theory of 
stratified sets. 

Let M be a smooth manifold, which we call the ambient manifold. Consider a 
singular subset V C M. Roughly speaking a stratification of V is a decomposition 
of V into a disjoint union of manifolds (strata) {V^} a such that strata of bigger 
dimension are attached to strata of smaller dimension in a "regular" way. 

"Regular" will obtain a precise meaning in a moment, but the most important 
property is that transversality to a smaller stratum implies transversality to an 
"attached" bigger stratum. Now we are going to describe the standard language of 
stratified manifolds and maps of stratified manifolds. This goes back to Whitney 



and Thorn fM, [Th 



Recall the Whitney Conditions (a) and (b). Condition (a) is similar to the notion 
of a p-stratification due to Thorn [|Th| d efined in the next subsection. We shall use 
ap-stratification to prove condition (|1C|). 

Consider a triple (Vp,V a ,x), where Vp, V a are C 1 manifolds, a; is a point in Vg 
and V/3 QV a \V a . 

Definition 28. A triple (V/3,V a , x) satisfies the Whitney (a) condition if for any 
sequence of points {xk} C V a converging to a point x £ Vg the sequence of tan- 
gent planes Tfc = T Xk V a converges in the corresponding Grassmanian manifold of 
dim V a -planes in TM and limTfc = r D T x Vp. 

Definition 29. A triple (Vfs,V a , x) satisfies the Whitney (b) condition if for any 
two sequences of points {xk} C V a , {yk} C V a converging to a point x £ Vp the 
sequence of "vectors" converges to a vector v e T X M which belongs to a 

limiting position of YimT Xk V a — T, i.e. v G r. 

Since condition (b ) is local one can think of M as Euclidean. This explains how 
to interpret the vector i^ZS^i • 
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It is easy to show that condition (b) implies condition (a). 

Definition 30. A locally closed subset V in the ambient manifold M is called a 
stratified manifold (set, variety) in M, if it is represented as a locally finite disjoint 
union of smooth submanifolds V a of M , called strata, of different dimensions in such 
a way that the closure of each stratum consists of itself and the union of some other 
strata of strictly smaller dimensions, and Condition (b ) of Whitney is satisfied. 

Any union of submanifolds satisfying condition of this definition 

(52) V = U a V a 

is called a stratification of V, and the submanifolds V a are called strata. A set V 
is stratifiable if there is a "nice" partition into strata. By a stratified manifold we 
mean a pair (V, V) consisting of a manifold V itself and a partition V = {Va}. 

4.6. Stratified maps and ap-stratification. Now we define a smooth map of a 
stratified manifold (V, V): 

Definition 31. Let (V, V) be a stratified manifold in an ambient manifold M, V C 
M , then a map f : V — > N is called C 2 -smooth if it can be extended to a C 2 smooth 
map of the ambient manifold M F : M — > N whose restriction to V coincides with 
f- 

A stratification V = U a V a stratifies a smooth map f : V — > M. k if the restriction 
of f to any stratum V a has constant rank, i.e., rank df\\r a (x) is independent of 

x e V a . 

A map G : L — > M is called transversal to a stratified set (V, V) ifG is transversal 
to each strata V a 6 V. 

By the Rank Theorem, if a stratification (V, V), V = {V a } a& i stratifies a smooth 
map P, then for each strata V a the number d a (P) = dim V a — rank dP\v a is well 
defined. 

Assume d a (P) > dp(P) for each Vp C V a \V a , i.e. nonempty level sets inside 
the bigger stratum V a have dimension d a (P) greater or equal to dimension of the 
level sets dp(P) in the smaller stratum Vp. We require that for any sequence of 
points {at} C P{V a ) converging to a point a 6 P(Vp), the nonempty level sets 
{P^ 1 (ak) nV a } approach the limiting level set {P~ 1 (a) O Va} "regularly" . In other 
words, we require that the level sets in the bigger stratum V a approach the limit 
level set in the smaller stratum Vp nicely. 

Definition 32. Let P : M — > N be a C 2 smooth map of manifolds, and let Vp and 

V a be submanifolds of M such that the restrictions P\v f3 to Vp and P\v a to V a have 
constant ranks R\rp(P) and Rv a {P), respectively. Let x be a point in Vp. 

We call the manifold V a ap -regular over Vp with respect to the map P at the 
point x if for any sequence of points {x n } C V a converging to x 6 Vp the sequence of 
tangent planes to the level sets = ker dP\y a (xk) converges in the corresponding 
Grassmanian manifold of (dim V a — Ry a {P))- dimensional planes to a plane t and 

(53) limfcer dP\v a (xk) —tD ker dP\vAx) 

Definition 33. AC 2 smooth map P : V — > TV of a stratifiable manifold V to a 
manifold N is called ap -stratifiable if there exist a stratification (V, V) such that the 
following conditions hold: 

a) (V, V) stratifies the map P (see definition 
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b) for all pairs Vp and V a from V such that Vp C V a \ V a the stratum V a is 
ap -regular over the stratum Vp with respect to P at point x for all x G Vp. 

The original definition of ap-stratification requires an appropriate stratification 



of the image also [Ma], but we do not require stratification of the image for our 
purposes. 

4.7. Relation between existence of ap-stratification and condition (|5|). In 

section 



1.5.1 we showed that the key to the proof of Theorem |26| is condition (10) 
(see Proposition Q). Now we are going to reduce the question whether condition 
( |To| ) is satisfied to the question whether an ap-stratification of the polynomial P 
exists. 

Let P = (Pi,P 2 ) : M w ^ I 2 be a nontrivial polynomial, V = P^ 1 ^) and 
Vb = (Pi,P2) _1 (0) be level sets. Assume that there exists a stratification (V, V) 
that stratifies the map P\v such that the zero level set Vb can be represented as a 
union of strata from V, i.e., Vq = U Qe / V Q . Denote this stratification of Vq by Vo- 
Recall that a map F : M. k — > is transversal to a stratification (Vb,Vb) if it is 
transversal to each strata V a G Vo- Associate to each level set V a , a^O a natural 
decomposition V a = {V a H V a } a ei- 

Proposition 2. With the above notation if a stratum V a G V \ Vo is ap-regular 
over a stratum Vp G Vo with respect to the polynomial P , then any C 2 smooth map 
F : R n — > R 2 transversal to (Vb, Vo) is also transversal to V a PI V a for any small a. 
This is equivalent to condition hit 



Proof Pick a point x in Vp C Vq and a point y G V a . Notice that ker dP\vp{x) is 
the tangent plane to the level set {P^ 1 (P(x))DVp} at the point x and ker dP\v a (y) 
is the tangent plane to the level set {P^ 1 (P(yj) n V a }. 

By condition ( p3| ) if a map F : X — > R N is transversal to ker dP\va(%) at a 
point x, then F is transversal to ker dP\v a (y) for any y G V a near x. 

Therefore, the condition "F is transversal to Vp at a point x" implies the condi- 
tion "P is transversal to V^ D V a for any small a" . This completes the proof. 

4.8. Existence of ap-stratification for polynomial maps. The existence of 
ap-stratifications is not a trivial question. There are some obvious obstacles. For 
example, let V C R™ be an algebraic variety and let P : R n — > R fc be a polynomial 
map. Assume that (V, V) stratifies P. If we have two strata V a and Vp so that 
V a lies "over" Vp C V Q \ V Q , then condition (^) can't be satisfied if dimension 
of the level sets d a (P) in the upper stratum V a is strictly less than that of dp (P) 
in the lower stratum Vp, i.e., dim/cer dP|v cv (y) < dim/cer dPlv^. In this case a 
plane ker dP|y (3 (x) of the lower stratum Vp should belong to a plane r of smaller 
dimension (see condition (|53|)), whi ch is impossible. Thorn constructed the first 



example when this happens [ GWPL ] . 
Thorn's example 

Consider the vector-polynomial P in the form P : (x,y) — > (x,xy). The line 
{x = 0} is the line of critical points of P. Outside of the line {x = 0} P is a dif- 
feomorphism. Therefore, the preimage of any point a ^ P _1 (a) is 0-dimensional. 
On the other hand, the preimage of is the line {x = 0}. 

Definition 34. Let us call an algebraic set V rank compatible with a polynomial P 
if there exists a stratification (V, V) which stratifies P and for any pair V a and Vp 
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from V such that Vp C V a \ V a dimensions of the levels dp (P) in the lower stratum 
Vp do not exceed dimensions of the level sets d a (P) in the upper stratum V a . 

It turns out that even if an algebraic set V is rank compatible with a polynomial 
P, then ap-stratification still does not always exist. Let us present an example 
with this property. The example below belongs to M.Grinberg. It seems that the 
existence of a counterexample was known before, but we did not find an appropriate 
reference. 

4.9. Nonexistence of ap-stratification. Let V — {(x,y,z,t) £ R 4 : x 2 — 
t 2 y + z} be the three dimensional algebraic variety and P : V — > R 2 be the natural 
projection to the last two coordinates, i.e. P : (x, y, z, t) — > (z, t). 

Proposition 3. With the above notations the set V is rank compatible with the 
polynomial map P and does not have ap-stratification. 

Proof Consider a rank stratification of V. Such a stratification consists of three 
stratum: V x = {x = t = z = 0}, V 2 = {t = 0, x 2 = z, x^ 0}, and V 3 = {t ^ 0}. 
On each stratum rankPlv; = i — 1. Level sets P~ 1 (t, z) — parabolas for t ^ and 
lines for t = 0. 

Show that for each point a = (0, a, 0,0) 6 V\ there exists a family of level sets 
such that at the point a the property ap-regularity of V3 over V\ fails. 

Consider the preimage of the curve {z = —at 2 } C R 2 . This is an algebraic 
variety of the form W a = {x 2 = t 2 (y — a)}. One can see that W a is the Whitney 
umbrella. The level x 2 = t^{y — a) is the parabola. As to — > this parabola tends 
to semiline x = t = z — 0, y > a. At the point a £ Vi the property ap-regularity 
of V3 over V\ clearly fails. This completes the proof of the Proposition. 

Let us mention a positive result on existence of ap-stratification. 



Theorem 35. [ Hirl | If V C R" is a semialgebraic variety and P : R" — > R is a 
polynomial function, then there exists an ap-stratification of (V,V) with respect to 
P. 



5. Existence of op-stratification. 

In this section we prove existence of ap-stratification in the special case we are 
interested in. As the Example || shows, the existence of a ap-stratification is a 
nontrivial question. In general, it does not exist. Unfortunately, the existence of a 
ap-stratification in our case does not follow from the classical results, so we need 
to prove it. 

Let R N and R fc be Eucledian spaces with the fixed coordinate systems x = 
(xi, . . . , xn) S R n and a = (ax, . . . , a^) £ R fc with N > k and a non-trivial vector- 
polynomial P : R N e R fe . Recall that P is a nontrivial if it has a point x <E R , 
where rank dP(x) — k. In what follows we call vector-polynomial by polynomial for 
brevity. 

Definition 36. Let m = (1, 7712, • ■ ■ , rrik) G 2^ and 8 > 0. We call the (m, 5)-cone 
Km,s the following set of points 

K m< s = {a= (at, . . . , a fc ) £ R fc : < a x < 5, 
! " 4) 0< |o i+ i| < |a 1 ...a j p+ 1 forj = l,...,k-l}. 
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Let m' = (1, m' 2 , ■ ■ ■ , m' k ) 6 ih. Define m' >- m // m' / m and > for all 
2 < i < k. We call the (m', S')-cone K m i j a refinement of the (m, 5)-cone K m j if 
m' >- m and S > 5' . 

Define the following sets 

(55) V m j,p = closure{P _1 (i ; s: mi 5)}, Vo, m ,p = ^s>oV m ,s,p 
Then one has 

Theorem 37. For any nontrivial polynomial P there exist an integer vector m £ 
Z" and positive 5 such that the following conditions hold 

a) the set Vq — Vo, m ,p (see $>o\)) is semialgebraic. 

b) the set V mt s.p consists of regular points of P, i.e. if b € V m .$_p, then the level 
set P _1 (6) is a manifold of codimension n. 

c) there exists a stratification ofVo by semialgebraic strata (Vq, Vq) satisfying the 
property: V m ^,p is ap-regular over any strata V a £ Vo with respect to P. 

In order to prove Theorem [37] we reformulate it in a convenient for us language. 
Let a e R fc . Denote by L a = P _1 ( a ) the level set of P. Recall that a e R fc is 
called a regular value if for any x 6 L a the rank of linearity of P is maximal, i.e. 
r&rikdP(x) = k. 

Definition 38. Let a,b e R k be values of P : R N e R k , B N C 1* be the unit 
ball centered at the origin, and 

do : B N x R fe G K, do(x, a) = inf \\x - y\\ 2 

y<£L a nB N 

(56) . n 

4:1 xle R, d (a, b) = sup d° {x, a). 

xeL b nB N 

Then the C° -distance between level sets L a D B N and L^ D B N 

D° P (a, b) = dist° c (L a n B N , L b n B N ) = - (d (a, b) + d (b, a)) . 

For any 1 < m < N denote by G m ' N the set of m-dimensional planes in the N- 
dimensional Euclidean space. G m ' N is so-called the grassmanian manifold. Below 
we introduce convenient for as distance in the grassmanian manifold G m,N . Now 
we define C 1 -distance between regular level sets L a and Lb in an appropriate for us 
way. Write P using coordinate functions P = (Pi, . . . , P k ) : R N -> R fe . If x e M N 
is a regular point of P, then gradients VPi(a;), . . . , VPfe(x) are linearly independent 
and span the space which is the orthogonal complement to the tangent space to the 
level set P^ 1 (P(x)). Define a Gramm-Schmidt orthogonalization operator: 

Definition 39. Let v\,...,Vk £ R fc be linear independent vectors. Define the 
Gramm-Schmidt linear operator by 

*:(v 1 ,...,v k ) = (vl,...,vl), where 

(57) * * {v2,vi) , V"^ ( v k, v j) 

j<k 



Remarks 0. The Gramm-Schmidt linear operator * has nothing in common 
with the asterisk operator used for the Khovanski reduction procedure in section 

1. Vectors {v±, . . . , v k } and {u*, . . . , v^} span the same k-dimensional space de- 
noted by L; 
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2. Vectors {y\, ... ,v£} form an orthonormal basis in the plane L; 

3. Let {Lt}t£(o,i] be a family of k-dimensional planes in M. N spanned by a 
family of vectors {vi(t), . . . , ^fc(i)}te(o,l] depending continuously on t. Consider 
{v*(t), . . . , Vfc(t)}t£(o,i] — ■ ■ ■ , Ufc(i)) as the family of orthonormal basis in 
{Lt} te (o i i]. Then sufficient condition that L t — > L in the grassmanian manifold 
G k,N is existence of an orthonormal basis {vl, . . . ,v^} of L such that 

58) . .\ 3 , ' V r -► as t -»• 0. 

Define the Gramm-Schmidt operator for the polynomial map P 

(59) {*(dP)i(x), *(dP) k (x)) = *(VPi(x), VP k (x)). 

Each vector *(dP)j(x) is given by the rational function in x. 

Let E C be the set of critical points of P. To measure C 1 -distance between 
two regular level sets we introduce the following function: Let x, y ^ E. Then 

R (r v] = v A «n-(4*(n-M) 2 \ 

(60) Pl £jv (*(nw>*(nw)(*(n-(y) ; *(n(i/))j 



(x,y) = \\x-y\\ +R P (x,y) 



Definition 40. Let a, b e R fc &e regular waZues ofP:R N ^> R fe , and L a = P -1 (a) 
and Lt, = P _1 (o) fee regular level sets. 

d° 1P :B N xM s \S^M, d? P (x,a) = inf Q P {x,y), 

y£L a nB N 

d hP : (R fc \ E) x (M' £ \ E) -► R, di,p(o,6)= sup d? P (x,a). 

x£L b nB N 

Then the C 1 -pseudodistance between regular level sets L a (~l P w and L;, n B is 
defined by 

(61) D l P {a, b) = dist c i (L a n B N ,L b nB N ) =^ (di, P (o, 6) + d ljP (6, a)) . 

Remark 3. FFe caH the function Dp(a,b) C 1 -pseudodistance, not C 1 -distance, be- 
cause it does not satisfy the triangle inequality. However, it satisfies the following 
triangle-like inequality 

(62) 2(Dp(a,b) + Dp(b,c)) > Dp(a,b). 

The reason we define ^-pseudodistance Dp(a,b) in such a way is because the 
function Dp(a,b) is algebraic (see Lemma^ below). 

The inequality ^6^) can be proven as follows. Let v,w € R w be vectors. Denote 
by Z(v,w) the angle between v and w. Direct calculation shows that 

k 

R P {x,y) =^sin 2 (Z(*(dP) j ( a; ),*(^)i(j/)))- 

3=1 

It is easy to check that 2 (sin 2 a + sin 2 /3) > sin 2 (a + (3) which is sufficient for the 
proof of the inequality ^62J . 

Now we can reformulate Theorem a7\ in the following way 
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Theorem 41. For any nontrivial polynomial P there exist an integer vector m = 
(1, 7712, • • ■ , TOfe) £ Z" and positive S such that the following conditions hold 

a) for any two values with the same first coordinate a = (i, a%, . . . , a^) and b = 
(t,b 2 , ...,bk) in K m . 5 

D° P (a,b) < t; 

b ) the same as in Theorem [g^ ; 

c) for any two values with the same first coordinate a = (t, a2, • • • , cife) and b = 
(t, 6 2 , ... , bk) in K m ,s 

Dp(a,b) < t; 

Let us show that parts a) and c) imply parts a) and c) of Theorem |37] respectively. 
Proof a) => from a) of Theorem 1371 Consider an algebraic curve of the form 
7(f) = (t J t"»2+i )i (m 2 +2)(m 3 +i) ) _ -) ^+2)...(m h +i))_ ne can check that 7(f) G 

K m .s for any t £ (0,(5). Denote by Vt,p = P _1 (7(i)). The set Uo<t<sVt,p is clearly 
semialgebraic set. By the Tarski-Seidenberg theorem the following set 

closure{U < t <,5V t ,p} \ {{J <t<5Vt,p} 

is semialgebraic. Since for any smooth curve j'(t) = (t, 72(f), ■ ■ ■ , 7fc(*)) S Km,8, t 6 
(0,(5) Hausdorff distance between the level sets Vt.p = -P -1 (7(f)) and VJ P — 
P _1 (7'(i)) is at most t, i.e. D p {^{t),^'(t)) < t. It implies that Hausdorff dis- 
tance between any two level sets of the form Vt.p and V/ P tends to as t — > 0. 
Therefore, 

closure{Uo<t<5Vt, P } \ {U 0<t < s Vt,p} = closurejp- 1 ^^)} \ {P^{K m j)}. 

This completes the proof of part a) . 

Proof c) from c) of Theorem |37]. Let us use notations of the proof of part 
a). By theorem ^ there is a stratification of Vo, m ,p such that the semialgebraic set 
{^o<t<sVt,p} is a ap-regular over Vo, m ,p. Indeed, let 7Ti : R fc — > K be the natural 
projection onto the first coordinate. Then a polynomial function p = 7ri o P : M. N — > 
R is well-defined and p _1 (f) = P _1 (7(i)). Application of theorem |35| to the map 

p : closure{U <t<5V" ti p} K 

gives existence of a required stratification. 

Since C 1 -distance between any two level sets of the form Vt.p — -P _1 (7(f)) and 
VI P = P^ij'it)) is at most t, i.e. D P (j(t),Y(t)) < t. It implies that (^-distance 
between any two level sets of the form Vt.p and V( P tends to as t — > 0. Therefore, 
ap-regularity of P~ l \K mt s) over Vb, m ,p follows from ap-regularity of {Uo<t<sVt,p} 
over Vb, m ,p. This completes the proof of part c). 

Before proving Theorem ^l] let us formulate a basic fact from elimination theory 

5.1. Elimination theory. Let C" 1 denote the m-dimensional complex space z = 
(z\, . . . , z m ) £ C m , m £ Z + . A set V in C m is called a closed algebraic set in C m 
if there is a finite set of polynomials F\ , . . . F s in Z\ , . . . , z m such that 

y(Fi,...,F s ) = {(z 1 ,...,z m ) £C m \ F j {z 1 ,...,z m ) = 0, l<j< s}. 

One can define a topology in C m , called the Zariski topology, whose closed sets are 
closed algebraic sets in C"\ This, indeed, defines a topology, because the set of 
closed algebraic sets is closed under a finite union and an arbitrary intersection. 
Sometimes, closed algebraic sets are also called Zariski closed sets. 
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Definition 42. A subset S of C m is called constructible if it is in the Boolean 
algebra generated by the closed algebraic sets; or equivalently if S is a disjoint 
union T\ U • •• U T/., where Tj is locally closed, i.e. Ti = T[ — T" , T- — a closed 
algebraic set and T" C T[ — a smaller closed algebraic. 

One of the main results of Elimination theory is the following 

Theorem 43. ( §jh\\ , ch.2.2) Let V C x C N be a constructible set and n : 
x C N -> O 1 be the natural projection. Then w(V) C is a constructible set. 

6. Proof of Theorem 

6.1. Existence of the (m, <5)-cone K m ,s of regular values of P (or Proof 
of Part a) of Theorem |4l| ). The set of critical values Ep C R fc of a nontrivial 
polynomial map P : R N — > R fc is an algebraic set of positive codimension. It follows 
from Sard's lemma for algebraic sets [Mi]. Suppose d : M. k — > R is a nonzero 
polynomial whose zero level set d _1 (0) D E. Fix coordinate systems in M N . By 
writing the linearization matrix dP : M. N — > R fc and considering (N — k + 1) different 
k x k minors one can calculate d explicitly. 

Lemma 2. For a nonzero polynomial d : R fc — ► R £/iere exists an integer vector 
m = (1, m2, . • • , mk) € Z?[ and <5 > smc/i iftai d does no< vanish on the (m. 5)-cone 

Remark 4. If Tip C d _1 (0) and Ep n A" m ,<5 7^ 0, i/iera i/iere exists x G A" mj( 5 suc/i 
</iai d(cc) = 0. T/iis shows that part a) of Theorem follows from this Lemma. 

Proof Let us prove the statement by induction in dimension fc. 
For k — 1 the level set d _1 (0) C R is a finite collection of points and Lemma is 
obvious. 

Without loss of generality assume d(xi, . . . , x k ) is not divisible by Xk- If d 
is divisible by Xh, then for some (3 G Z + one can decompose d(x\, . . . , x k ) — 
Xud(xi, . . . , Xk) so that d(:ri, . . . , Xk-i, 0) is not identically zero. If for some m € l\ 
and 6 > the (m, 5)-cone K m s does not intersect zero locus d~ 1 (0), then K m s 
does not intersect zero locus d _1 (0) = d _1 (0) U {xt — 0} too. 

With the assumption of indivisibility by x k the following set E^T 1 = d _1 (0) n 
{Xf. = 0} C R fe_1 is of a positive codimension. By inductive hypothesis there exist 
an integer vector nih-i £ Zjj. -1 and 5i--i > such that the (m.fc_i, <5fc_i)-cone 
if' 1-1 r C R fc_1 has empty intersection with Ep -1 . 

Let a = (ai, . . . , a^) £ and |a| = • a^. Write d(x) = X) Q ez fc a a x a Denote 



by degd = max{|a| : a G ih, a a ^ 0}. Put m k = degd+ 1. 



r. 4 nl p- 

Proposition 4. With the above notations there exists 8 > such that for m = 
(m.fc_i,TOfe) 6 Z^_ i/ie (m, 8)-cone K m s does not intersect zero locus d _1 (0). 

Proof Put d(0) = otherwise the proposition is trivial. Write 

dcg d 

d(xi, ...,X k )= ^ dj(xi, . . .,x k -!)x 3 k . 

3=0 

Without loss of genericity one can assume that ao(xi, . . . , x^-i) does not van- 
ish on the (m^_i, <5fc_i)-cone K^ 1 s . If not, then apply Lemma || and re- 
fine 1 $ k 1 to a required size. By the definition of the (m, o")-cone K m g 



3(5 



Bifurcation of polycyclcs 



the condition (xi, . . . , Xfe) G K m $ implies that (xi, . . . , Xk-i 

) e K ml 5 fc _ and 

< x fe < (xi . . . Xfc_i) mfc . Put x k = A(xi . . . Xfe_i) mt: with A G (0, 1). It is easy to 
check that 

(63) d(xi, ...,x k ) = a (xi, . . .,Xk-i){l + p(xi, . . -,Xk-i, A)), 

where p is such a polynomial that p(0, A) = 0. Indeed, the choice of rrifc is 
such that (x\ . . . Xk-i)" lk = ao(xx, . . . , Xk-i)q(xi, . . . , Xk-x) for some polynomial 
q(xi, . . . , Since p(0, A) = for a sufficiently small 8, any (xi, . . . , Xk-x) G 

^nT_ 5 _i' and an Y G [0,1] the following inequality holds |p(xi, . . . , x k -i, A)| < 
1/2. This shows that d does not vanish on the (m, <5)-cone and completes the 

proof of the Proposition. 

As we pointed out above the Proposition implies Lemma |^. 

6.2. Reduction to an optimization problem (or Proof of parts a) and c) of 
Theorem |37|). Let P — (P%, . . . , Pf.) : M. N — > R fc be a nontrivial polynomial with 
N > k given by its coordinate functions and an (m, 8) — cone -K7 m ,5) C Im P(R N ) C 
R fc be a cone of regular values of P. Existence of such a cone is proven in the 
previous section. Recall that E C R N denotes the set of critical points of P and 
Qp(x, y) is defined in (|6^). The function Qp(x, y) is a rational function symmetric 
with respect to permutation of x and y. It is defined to measure C^-distance 
between level sets (see remarks after definition [39]) . The singular set of Qp belongs 
to (E x R N ) U (R w x E). Recall that P w = {x : a;? < 1} C M^. Introduce 

functions r(x) = 1 — ^j =1 xf . 

Assume that the restriction of P to the boundary S N = dB N = {x : ||x|| = 1} 
has only the regular values in the (m, 5)— cone K^ m ^y Indeed, regularity of P|s« : 
S N — > R* is equivalent to regularity of the polynomial map (P, r) : — » R fe x R 
given by (P,r)(x) = (P^^ili^i)- Existence of an (m, 6)— cone of regular values 
of the map (P, r) follows from Lemma ^. 

Lemma 3. With the notations above let a T1 , a T2 G if ( m ,5) & e points with the 
same first (fc— 1) coordinates, i.e. a Tl — (a fe ~ 1 ,ri) and a T2 — (a , T2). Then there 
exists a polynomial R(a k ^ 1 , Ty, T2, c) m variables and 
c G R smc/i i/iai 

(64) i?(a fe - 1 ,ri,r 2 ,rfi :P (a Tl ,a T2 )) = 0. 

Moreover, R{a k ~ 1 , r, r, 0) = 0. 

Proof Recall that {*(<iP) ;) -(x)}|_ 1 form an orthogonal basis in the orthogonal 
complement to the tangent plane to the level set P _1 (P(x)) at the point x (see 
(59)). Let us make several remarks about the rational function i?p(x, y) defined by 
©• 

1 . If a T = (a k 1 , r) e K m< s is a regular value for the map (Pi , . . . , P& ) : R N — ► R fc 
for some r, then a k ~ 1 G R*^ 1 is a regular value for the map (Pi, . . . , Pfc-i) : R^ — » 
R fc_1 ; 

2. If a' 1 '- 1 G R fe_1 is a regular value for the map (Pi, ... , P k -i) ■ B N -> R*^ 1 , 
then there is a positive e(a fc_1 ) > such that for each point x G (Pi, . . . , P/ C _i)~ 1 (a' £-1 ) 
and each 1 < j < k — 1 

(65) HdPUx), *(dP)j(x)) > 6(a fe - 1 ) > 0. 

This follows from compactness of B and regularity of the value a ; 
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3. Since we consider only those r that a T belongs to the (m, i5)-cone K m s of 
regular values of P there exists a positive constant e(a k ~ 1 , r) > such that for each 

x g P -1 (or) 

(66) (*(dP) k (x),*(dP) k (x)) > e(a k -\r). 

This shows that Q(x, y) restricted to P^ 1 {K m ,s) x P~ 1 {K m: $) is a smooth function 
of x and y. 

Consider irreducible representation of the rational function Q(x,y) as a ration 
of two polynomials Q(x,y) = ■ Because of remarks 2 and 3 S(x,y) ^ for 

each pair (sc,y) £ P -1 (.Km,*) x P'H-^m.s)- 

Now notice that we deal with smooth objects: smooth level sets P~ 1 {a T ) and 
the smooth function Qp(x,y). Notice that d\ P (x 1 a T2 ) is an extremal value of 
the function Qp(x 7 y) provided that P(y) = a T2 . Similarly, di } p(a Tl , a T2 ) is an 
extremal value of the function d\ P (x,a T2 ) provided that P{y) — a Tl . To find 
all extremal values of a smooth function on a smooth manifold one can use the 
Lagrange multipliers method. We prove that functions d\^p and d® P are algebraic 
functions. 

The key point of the Lagrange multipliers method is that at an extremal point of 
Qp(x, y) under the condition P(y) = a T2 the gradient V y Q P (x, y) can be expressed 
as a linear combination of gradients VPi(y), . . . , VPfc(y), and Vr(y). The gradient 
of Qp{x,y) has the form 

V y Qp(x, y) = V y ( ^y~^T ) = ( S p( x > y)^y T p( x , v) ~ t p(x, y)V y S P (x, y)) S P 2 (x, y). 



Since S\p-ifK m s )xP- 1 (K m e ) / " we can rewrite the Lagrange system in the 
following form 

' Sp(x,y)VyT P (x,y) - T P (x,y)V y S P (x,y) + 



(67) 



0, 



+Sp(x) [E-=i AiVPjfc) - X k+1 Wr(y) 
P{y) - a T2 

T P (x, y) - cS P (x 1 y) = 0, 
K A fe+1 r(y)=0. 

Important that all equations are polynomial and we can apply elimination theory! 
Notice that the last equation is responsible for an extremal point y which might 
belong to the boundary dB N . If a critical value belongs to the boundary, then 
r(y) = and \ k+ i\/r{y) is not zero and the gradient V y Qp{x,y) should be ex- 
pressed as a linear combination of k + 1 vectors VPi(y), . . . , \7P k (y), and Vr(y). If 
a critical value does not belong to the boundary, i.e. r(y) ^ 0, then Afe+i = and 
A fc+ iVr(y) = 0. 

Complexify the system ([37]), i.e. consider the system (|6T|) for 

(x,y,X,a k - 1 ,T 2 ,c) G C N x C N x C^ 1 x C*" 1 x C x C. 

It defines a constructible set, denoted by V, in £ 2N + 2k + 2 _ L e t us eliminate vari- 
ables {yi]f=i by projecting V along the corresponding (JV + k + 1)- 
dimensional (A, y)-plane. The result of projection is a constructible set W in the 
space (x, a k ~ 1 , T2,c) g C Ar+fe+1 . By the construction a point (x, a k ~ 1 , t<z, c) belongs 
to if some value of y the following conditions hold: P(y) — a T2 , Qp(x,y) = c, 



and y is the critical point of Qp restricted to P 1 (a 



T 2 , 



:!8 
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The constructible set W has dimension N + k. Indeed, consider a polynomial 
function p x : P~ 1 (a T2 ) — > R dehned by p x (y) — Qp{x,y). By Sard's lemma for 



algebraic sets Mu| critical values of p x form an algebraic set of positive codimension 
in M. Therefore, the set of critical values consists of a finite number of points and 
a finite number of possible c so that (x, a k ~ 1 , r%, c) G W. Thus, dimW equals 
dimension of (x, a fe " 1 , T2, )-plane, i.e. dim W = N + k. 

Since W is constructible and has codimension 1 there is a non zero polynomial 
R(x,a T21 c) such that W C R~ l (0). By the definition ( |6l| ) of d\ P {x 1 a T2 ) and by 
the construction 

(68) a Tl , a T2 , d° p(x, a T2 )) = 0. 

In order to prove that the function di ! p(a T1 , a T2 ) defined by ( |6l| ) is also alge- 
braic, calculate critical values of d\ P (x, a T2 ), provided P{x) — a Tl . By the implicit 
function theorem the gradient V^dJ p(x, a T2 ) can be expressed in terms of partial 
derivatives of R(x, a T2 , c) by the following way 

(69) d Xj d° 1P (x, a) = d c R(x, a, c) (d Xj R(x, a, c)) 

for (a, c) = (a T2 , d® p (x, a T2 )) an d provided that d Xj R{x, a T2 , d\ P {x, a T2 )) ^ for 
all 1 < j < N. Fix a = a T2 and consider x outside of the union of algebraic sets B = 
ujLi{x : d Xj R(x,a,c) = 0}. Then d° 1P { x, a) is a smooth function in x. Application 
of the Lagrange multipliers method shows that at an extremal point of the function 
di P (x,a), provided P(x) = a Tl , the gradient V^d" P (x, a) can be represented as 
a linear combination V x d® P (x,a) — ^ - =1 Xj VPj(x) — Xk+iVr(x). Plugging in 
the expression for V x d^ p (x.a) in terms of 9 c i?(a;, a T2 , c) and d Xj R(x,a T2 ,c) for 
j = 1, . . . , N we can present a Lagrange multiplier system in the following form 

d Xj R(x, a T2 , c)d c R(x, a T2 , c) = 

= EjLi Aj VPjOx) - A fe+ iVri(a:) (d Xj R(x, a T2 , cfj 

( 70 ) { R(x,a T2 ,c) = 0, 

P(x) = Or, , 
A fe+ ir(a;) = 0. 

Again the system ( f70|) consists of only polynomial equations and we can apply 
elimination theory. Consider this system for 

{x,X,a k - l ,T U T 2 ,c) £ C N x C fc+1 x C^ 1 x C x C x C. 

It defines a constructible set, denoted by Vi, in C Ar+2fc+3 . Let us eliminate vari- 
ables {A 3 -}£tj, {^il^j by projecting V along the corresponding (N + k + 1)- 
dimensional (A, x)-plane. The result of projection is a constructible set W± in 
the space (a k ~ 1 , t\, t 2 , c) e C A ' +2 . 

Similarly to the arguments for the constructible set W one can show that Wi has 
dimension k. Since W\ is constructible and has codimension 1 there is a nonzero 
polynomial R(a k ~ 1 ,ti,T2,c) such that W\ C i? _1 (0). By the definition ( |6l| ) of 
di,p(a Tl , a T2 ) and by the construction 

(71) R(a Tl ,a T2 ,di,p(a Tl ,a T2 )) = 0. 
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By the construction if a Tl = a T2 , then R(a Tl , a Tl , 0) = 0, because in this case 
both level sets are the same and C 1 -distance between them must equal zero. This 
completes the proof of Lemma [| 

Lemma 4. With the notations above there exists a refinement (m', S')-cone K m > $> C 
-Km, £ such that for any pair of points a Tl = (a ^ti) and a T2 = (a ,T2) from 
K m ',S' 

(72) D 1 P {a Tl ,a T2 )<{a 1 ...aik~l) 2 , 
where a fc_1 = (ai, .. . , a fc _i) G R fc_1 . 

Proof It follows from Lemma |^ that there is a polynomial R(a Tl ,a T2 ,c) such 
that R{a Tl , a T2 , di t p(a Tl , a T2 )) = and R{a, t, t, 0) = 0. 

For a Tl , a T2 belonging to the (m, <5)-cone K m g of regular values of P the function 
di,p{a Tl , a T2 ) depends continuously on t\ and t 2 . Let us rewrite R{a Tl ,a T2 ,c) in 
the form R(a k ^ 1 , Ti, T2, c). Recall that in our notations a T = (a^ -1 ,? - ). 

Suppose for determiness that t\ > Notice that each sufficiently small positive 
root Cj(a k ~ 1 , Ti, r 2 ) is increasing in n and decreasing in t 2 in a neighborhood of 
(a k ~ 1 ,T 1 ,r 2 ) =0. Therefore, c J (a fe_1 , n, 0) > c-,(a fe ~ 1 , n, r 2 ) 

Denote R(a k ~ 1 , t, c) = R(a k ^ 1 , t, 0, c). Let us show that for some sufficiently 
large positive integers ml and m" if < r < (a^ 1 . . . a k z\c) m and < c < 
(a* -1 . . . a k k z\) m '' the following decomposition holds: Put c = ^(a^ 1 . . . a k k z\) m '' 
and r = A(a^ 1 . . . a k k z\c) m ' with p, A e (0, 1). Then 

(73) R{a k -\ t, c) = roCa^Xl + Ql (a*" 1 , p))(l + g 2 (a"" 1 , p, A)), 

where ro, ?i, 92 are polynomials in their variables. Indeed, apply the same argu- 
ments as we used to prove ( |63] ) to the polynomial 

dcgft 

R(a k - 1 ,T,c)= J2 Rj{a k ~ 1 ,c)r j + flo(a* -1 ,c). 
3=1 

Then apply the same arguments to 

dog Ro 

R (a k -\c)= r M k ^) c3 +r (a k ~\c). 

3 = 1 

Notice that i?(a fc-1 ,0, 0) = implies that <?i(0,p) = g 2 (0,p, A) = 0. Therefore, for 
a sufficiently small 6' > and any (a k ~ 1 ,5 l ) 6 K m ^' polynomials qi(a k ~ 1 , p) and 
q2(a k ~ 1 ,p, A) are sufficiently small and R(a k ~ 1 ,r, c) equals if and only if r (a k ~ 1 ) 
equals 0. 

By Lemma cone there is a refinement (m',5')-cone K m i s> C K m .$i such that 
ro(a k ~ 1 ) does not vanish on K m r 

Now put m' k — m'm" . As we have just shown all sufficiently small positive roots 
c j (a k - 1 ,T,0) < t- 1 /™'^- 1 . . .a k z\) provided that 

— l/m'/k — 1 _fc— 1\ ^ („k—X „k—l\m" 

T> {a 1 ...a k _ 1 )<(a 1 ■■■a k _ 1 ) . 

This condition is satisfied for any < r < (a^ -1 . . . a k z{) m ' m " ■ This shows that all 
sufficiently small positive roots 

Cj {a k -\ r, 0) < (a^ 1 . . . a k k z\) m " +1 < {af 1 ■ ■ ■ 4z\) 2 . 
This completes the proof of the Lemma. 
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Let us complete the proof of part c) Theorem ^l] by the following inductive 
arguments. 

Consider a sequence of positive integers mi, . . . , mfe. Let m = (1, mi, . . . , m^) 
and 5 > 0. Define a sequence of polynomials associated to this sequence, defined 
by their coordinate functions: 

P? = Pj-(Pi...Pj-i) m >, J = 2,...,fc 

(74) 

P s = (P°, . . . ,P?,P S+ 1, ■ ■ ■ ,Pk-s), s = 2, . . . ,k 

Define the restriction of the (m. <5)-cone K m j to the s-dimensional plane, denoted 
by s , generated by the first fc-coordinates by the following way: 

K s m S = {a s = (oi, . . . , a s ) G R s : < a x < 5, 

< |a i+i | < \ai . ..a \ m i +1 for j = 1, . . . , s - 1}. 

It is shown above that there is an (m, <5)-cone K m g such that any point (0, a k ~ 1 ) G 
R x K^ s is a regular point for the polynomial P k ~ 1 . Therefore, one can apply 
Lemmas |[ |3[ and ^ and show that there is refinement of K^g, denoted the 
same, such that for any two points a^ 1 = (0, a fc_2 ,Ti) and a^ _1 (0, at-2, T2) from 



x^ 1 



(76) D 1 ^- 1 ,^- 1 ) < ( ai ... afe _ 2 ) 2 . 

By induction one can show that there is a refinement an (m, <5)-cone K m s such 
that for any two points a£~ s = (0, a fc_s_1 , t\) and a^~ s (0, ak-s-i, T2) from the 
restriction cone Kz~? such that 

(77) D 1 P s(4- s ,a k r - s )<(a 1 ...a k . s ^ 1 f. 

Notice that for any 1 < s < k level sets of the polynomial P s correspond to level sets 
of the initial polynomial P. Combining this with all estimates for D 1 pa (a k .~ s , a k ~ s ) 
and the triangle-like inequality (p2|) one can show that part c) of Theorem [ll] holds 
true. Part a) of Theorem [ll] follows from part c) because Qp(x,y) > \\x — y\\ 2 , 
which implies that Dp(a, b) > Dp (a, b) for any pair a, b G K mt $. 
This completes the proof of Theorem 
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